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Abstract 


The  ATS-6  radio  beacon  experiment  carried  out  at 
Aberystwyth  (52.42  N,  4.05°W)  is  described.  A preliminary 
analysis  of  the  data  is  described  concentrating  on  hourly 
values  of  total  electron  content  (N.j,) , Faraday  content  (Np) , 
protonospheric  content  (Kp)  and  layer  shape  F-factor.  The 
monthly  median  behaviour  of  the  parameters  is  discussed 
together  with  the  effects  of  observational  and  analysis 
errors  on  the  accuracy  of  the  protonospheric  content  values. 

The  correction  for  layer  height  changes  has  been  found  to  be 
of  much  less  importance  for  European  stations  observing 
during  Phase  II  ATS-6  operations  chan  for  measurements  made 
in  America  during  Phase  I. 

The  dynamic  interchange  of  plasma  between  ionosphere 
and  protonosphere  has  been  studied  by  estimating  hourly  values 
of  average  integrated  net  fluxes  from  the  temporal  gradient 
of  protonospheric  content.  The  protonospheric  flux  data  have 
been  interpreted  qualitatively  in  terms  of  diffusive  interaction 
between  ionosphere  and  protonosphere  in  both  local  and  conjugate 
hemispheres.  Added  confirmation  is  obtained  by  consideration 
of  contrasting  protonospheric  data  from  an  American  sector 
station  which  can  also  be  related  to  ionospheric  behaviour  in 
the  two  hemispheres. 

The  Aberystwyth  observations  indicate  that  the  inner  flux 
tubes  intersected  by  the  slant  path  are  filled  in  a matter  of 
hours  almost  every  day  and  the  average  behaviour  of  Np  reflects 
a diurnal  partial  draining  but  complete  refilling  of  tubes  up 
to  L 't'  4.  Case  studies  show  that  following  a geomagnetic  storm 
depletion  some  7 days  are  required  to  saturate  these  tubes, 
but  approximately  14  quiet  days  are  necessary  to  reach  absolute 
dynamic  equilibrium  along  the  slant  path. 
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1.  INTRODUCTION 


The  ATS-6  satellite  was  launched  into  geostationary  orbit  in 
May  197A  and  was  stationed  at  longitude  9A°W  until  mid-1975  when 
it  was  moved  to  35°E  becoming  available  to  observers  in  Europe. 
Among  other  experiments  the  satellite  carried  a multifrequency 
radio  beacon  designed  for  several  ionospheric  investigations.  Of 
specific  interest  to  the  present  study  were  transmissions  which 
enabled  absolute  measurements  to  be  made  of  electron  content,  a 
parameter  of  importance  in  assessing  the  ionospheric  delay 
correction  in  very  high  frequency  communication  and  navigation 
systems.  Two  types  of  measurement  possible  with  the  ATS-6  radio 
beacon  are  of  concern  here.  In  the  first  the  total  electron 


content  (N^)  between  satellite  and  ground  can  be  determined  by 

comparison  of  the  phase  of  a IMHz  modulation  on  a VHF  carrier  with 

that  of  an  identical  modulation  on  a coherent  VHF  carrier.  Secondly, 

a weighted  estimate  of  the  ionospheric  electron  content,  the  so- 

called  Faraday  electron  content  (N  ) , restricted  to  a measure  of 

r 

the  integrated  plasma  below  some  arbitrary  height  at  the  top  of 
the  ionosphere,  can  be  obtained  from  the  rotation  of  the  polari- 
sation plane  of  the  VHF  carrier  transmission.  Absolute  deter- 
mination of  the  total  and  Faraday  electron  contents  to  a high 
accuracy  enables  the  so-called  protonospheric  electron  content 
(Np)  along  the  slant  path  to  be  determined  from  the  difference. 

The  present  work  is  concerned  with  preliminary  results  from 
ATS-6  observations  made  at  Aberystwyth  (52.42°N,  4.05°W)  from 
November  1975  to  July  1976  when  the  satellite  was  situated  at  35°E. 

As  particular  emphasis  is  placed  in  this  report  on  aspects  of  the 


li'otono.spheric  electron  content  measurement  the  description  of  the 
experiment  and  the  results  obtained  are  preceded  by  a short  review 
of  the  literature  relevant  to  plasmispheric  processes  and  the 
dynamic  coupling  between  the  ionospheric  F2-region  and  the  overlying 
protonosphere . 

2.  LITERATURE  REVIEW 

In  discussing  the  complex  interaction  between  ionosphere  and 
p’otonosphere  several  factors  should  be  kept  in  mind.  For  example, 
the  plasma  is  constrained  to  diffuse  along  the  geomagnetic  field 
lines,  the  configuration  of  which  on  a global  scale  is  complicated 
and  in  practice  subject  to  ill-defined  changes  in  both  space  and 
time  (HESS,  1968).  Secondly,  the  ionosphere  and  protonosphere  or 
exosphere  are  dominated  by  different  ions  0^  and  respectively 
with  a diffusive  barrier  between  them  (HANSON  and  ORTENBERGER,  1961). 
Thirdly,  the  effective  vertical  motion  of  the  ionospheric  plasma 
is  controlled  by  thermospheric  winds  and  electric  fields  whose 
patterns  on  a global  scale  are  complex.  It  is  thus  perhaps  not 
surprising  that  despite  considerable  interest  in  the  problem  for 
some  time  relatively  few  direct  measurements  have  been  made  of 
the  relationship  between  ionosphere  and  protonosphere  and  reports 
have  tended  to  concentrate  on  individual  events  with  the  data  too 
sparse  for  statist.’cal  study. 

Observations  of  whistlers  were  used  in  early  studies  to 
estimate  geomagnetic  flux  tube  content  and  the  associated  plasma 
fluxes  between  ionosphere  and  protonosphere.  PARK  (1970)  used 
whistler  measurements  of  tube  content,  for  a period  of  several 
days  after  a moderate  magnetic  storm,  and  from  the  temporal 
changes  estimated  the  magnitude  of  the  upwards  flux  of  electrons 
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from  the  ionosphere  during  daytime.  CHAPl’ld.L  et  al . (1971)  obtained 
confirmation  of  this  flux  rate  from  in  situ  measurements  of  density 
by  means  of  a light  ion  mass  spectrometer  aboard  the  OGO-5  satellite. 

An  early  attempt  to  obtain  continuous  measurements  of  protono- 
spheric  electron  content  along  the  line  of  sight  to  a geostationary 
satellite  was  reported  by  ALMEIDA  (1970).  Using  transmissions  from 
the  ATS-3  satellite  the  Faraday  content  and  the  total  content  were 
estimated,  the  latter  by  the  differential  phase  path  technique, 
allowing  the  protonospheric  contribution  to  be  obtained.  The 
results,  discussed  more  fully  by  ALMEIDA  (1973  and  74),  show  that 
while  in  principle  the  protonospheric  content  can  be  measured, 
apparently  anomalous  results  are  on  occasions  obtained  with 
minimising  during  daytime  when  the  upwards  flux  of  electrons  from 
the  sunlit  ionosphere  would  be  expected  to  be  greatest.  Similar 
diurnal  variations  were  obtained  from  a restricted  data  set  by 
KLOBUCHAR  and  KIDD  (1972),  for  ATS-3  observations  again  in  the 
Western  Hemisphere.  These  workers  suggest  that  the  anomaly  may 
be  an  observational  effect  arising  from  the  choice  of  a constant 
height  for  the  estimate  of  the  so-called  M-factor  used  in  the 
Faraday  content  estimation,  a point  discussed  more  fully  by 
ALMEIDA  (1974).  Semi-empirical  model  studies  by  SMITH  (1970) 
and  independently  by  KERSLEY  and  SAMBROOK  (1971)  and  SAMBROOK  (1974), 
combining  geostationary  satellite  Faraday  rotation  measurements 
with  incoherent  scatter  data  supplemented  by  a model,  implied 
that  the  protonospheric  contribution  to  the  total  content  could 
be  typically  about  10%  by  day  rising  to  perhaps  40%  at  night. 

The  idea  of  the  plasmasphere  acting  as  a reservoir  for  ions 
which  can  flow  to  the  ionosphere  when  the  F-region  pressure  is 
low  and  vice-versa  has  been  discussed  by  BANKS  and  DOUPNIK  (1974). 
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The  plasma  transport  between  ionosphere  and  protonosphere  implies 
an  ion  transfer  between  tlic  conjugate  hemispheres.  MAYR  et  al.(1972) 
using  model  calculations  obtained  an  estimate  of  tlic  flux  between 
the  hemispheres  to  account  for  aspects  of  ionospheric  behaviour. 

Theoretical  estimates  of  the  ionospheric/protonospheric  flux 
have  also  been  carried  out  by  solving  the  continuity  equation, 
using  composition  data  from  the  incoherent  scatter  technique 
(HO  and  MOORCROFT,  1975).  The  density  and  vertical  velocity  of 
the  different  plasma  constituents  measured  by  this  method  have 
been  used  to  estimate  plasma  flux  (EVANS,  1975  a and  b) . This 
method  is  perhaps  the  most  direct  way  of  estimating  the  vertical 
ionisation  flux  but  it  suffers  from  the  usual  operational  limitations 
to  only  a few  day  observations  per  month, while  in  addition  Che 
transition  height  from  0^  to  H is  often  much  higher  than  the 
maximum  height  at  which  the  incoherent  scatter  technique  produces 
a reasonable  signal  to  noise  ratio. 

The  multifrequency  radio  beacon  aboard  ATS-6  has  made  available 
a new  technique  for  the  determination  of  the  protonospheric  electron 
content  along  the  path  to  a geostationary  satellite.  The  overall 
description  of  the  ATS-6  radio  beacon  experiment  and  associated 
ionospheric  observations  have  been  discussed  by  HARGREAVES  (1970) 
and  DAVIES  et  al.  (1972).  Estimates  of  the  sensitivity  of  the 
technique  to  different  ionospheric  and  protonospheric  electron 
density  distributions  were  made  using  model  studies  by  HARGREAVES 
and  HOLMAN  (1972). 

Some  early  results  of  ATS-6  observations  at  Boulder,  U.S.A. 
were  reported  by  DAVIES  et  al.  (1975  a and  b)  while  a complete 
set  of  data  for  this  station  covering  the  period  July  1974  to 
May  1975  has  been  published  by  FRITZ  (1976).  DAVIES  et  al.  (1976) 
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have  linked  this  data  ?.o  model  calculations,  but  a feature  of  the 
results  wh'ch  passes  largely  without  comment  is  a diurnal  variation 


of  the  monthly  median  protonospheric  content  v.'hich  maximises  at 

night  during  the  winter  months.  SOICHER  (1975  and  1976a)  reports 

ATS-6  measurements  of  protonospheric  content  for  periods  of  a few 

days  one  of  which  shows  enhanced  values  followed  by  an  extended 

depletion  following  a sudden  commencement  storm.  SOICHER  (1976b) 

contrasts  the  behaviour  of  protonospheric  content  for  ATS-6 

observations  in  U.S.A,  South  America  and  Northern  Europe  while 

SOICHER  (1976c)  presents  multi-plots  of  data  for  Fort  Monmouth,  U.S.A 

for  four  months  in  1975.  ATS-6  measurements  of  N have  been 

P 

compared  to  model  simulations  by  POLETTI-LIUZZI  et  al.  (1976)  for 
two  periods  in  1976,  one  a long  string  of  quiet  days  and  the  other 
including  two  major  geomagnetic  storms.  The  diurnal  variations 
of  protonospheric  content  for  this  American  sector  station  are 
small  and  show  some  evidence  for  a daytime  minimum,  a result  again 
passed  over  without  comment.  By  contrast,  HARTMANN  et  al.  (1976), 
with  observations  from  Europe  in  a paper  primarily  concerned  with 
travelling  ionospheric  disturbances,  show  a daytime  maximum  in 
protonospheric  content. 

The  present  report  concerns  ATS-6  observations  from  Europe 
where  the  monthly  median  behaviour  of  the  ionospheric  and  protono- 
spheric contents  is  discussed  and  compared  to  observations  from  the 
Western  Hemisphere  and  consideration  given  to  the  plasma  dynamics 
contributing  to  the  observed  variations.  In  addition,  detailed 
attention  is  given  to  certain  individual  periods  of  quiet  and  enhanced 
geomagnetic  activity. 
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3.  THEORY 

a.  The  Basic  F^quacions 

Consider  a linearly  polarised  electromagnetic  wave  of  frequency 

f,  characterised  by  two  circularly  polarised  components,  tlie 

ordinary  (o)  and  extraordinary  (x)  modes,  propagating  along  a small 

path  increment  ds  in  an  ionised  medium  in  the  presence  of  a 

magnetic  field.  The  refractive  indices  will  be  given  by 

u = c/fX  so  that  the  phase  changes  of  the  two  modes  deb  and 
o,x  o,x  o 

d<))^  will  be  given  in  the  usual  notation  by 


2ir  w , 

d<t  = T — ds  = — y ds 

o.x  A c o,x 

o,x 


1. 


The  refractive  indices  are  obtained  from  the  Appleton  Hartrec 
equation,  which  for  QL  propagation  in  a collisionless  medium, 
approximations  relevant  to  the  present  work,  reduces  to 


^o,x  " ^ " 2(1±Y  1 


where  X = 


Ne= 


eB. 


■y  and  Y,  = — — 
E mw  L mw 


in  the  usual  notation  with  N 


the  free  electron  density,  the  component  of  the  magnetic  field 
in  the  propagation  direction  and  f the  longitudinal  electron 

Li 

gyrof requency . 

For  propagation  along  a path  s the  polarisation  rotation 
of  the  linearly  polarised  w?  v^e  will  be  given  by  half  the  accumulated 
phase  difference  between  the  modes,  that  is 


UJ 

n = J (y  - y Jds 

2c  s o X 


so  that  by  substitution 


n = - / XY,  ds 

c s L 


that  is 
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where  k has  a numerical  value  4.840  x 10  for  £1  in  degrees, 
with  the  frequencies  in  MHz  and  distances  in  metres. 

Considering  now  only  the  ordinary  mode  circularly  polarised 
component,  its  accumulated  phase  change  along  the  path  s will 
be  given  by 


which  by  substitution  becomes 


Lt 


that  is 


— s - k ; v -V  ds 
c s (f+fj^) 


The  first  terra  on  the  right  hand  side  of  this  equation 
represents  the  geometrical  phase  path  while  the  second  term 
corresponds  to  the  phase  delay  resulting  from  the  ionosphere. 
That  is  the  ionospheric  phase  delay  at  frequency  f is  given  by 


b.  Application  to  the  ATS-6  Experiment 


(i)  Faraday  content 

The  Aberystwyth  receiving  system  measures  the  phase  difference 
between  the  ordinary  and  extraordinary  circular  modes  of  the 
140.056MHz  carrier  allowing  the  Faraday  rotation  angle  (0)  to  be 
determined.  If  a mean  value  of  f.  can  be  determined  representative 

Xj 

of  conditions  along  the  ionospheric  path  then  Equation  3 can  be 


written 


where  N^  is  the  so-called  Faraday  content.  Because  of  the 

weighting  of  the  geomagnetic  field  N can  be  considered  as  a 

r 

measure  of  the  ionospheric  electron  content  up  to  some  arbitrary 
height  usually  between  2000  and  3000  km  (SAMBROOK,  1974)  . 

To  estimate  a mean  value  of  f^^  it  is  necessary  to  know  how 
this  parameter  varies  along  the  path'.  The  computation  has  been 
carried  out  using  a spherical  harmonic  model  of  the  geomagnetic 
field  with  coefficients  based  on  CAIN  and  SWEENEY  (1970)  but 
up-dated  to  epoch  1974  (DAVIES  1975).  Table  1 displays  an 
output  from  the  program  listing  values  of  several  geometrical 
and  geomagnetic  parameters  relevant  to  the  ATS-6  to  Aberystwyth 
geometry  at  50  km  vertical  intervals  through  the  ionosphere  to 
a height  of  2000  km.  The  values  of  f^^  over  this  height  range 
are  plotted  in  Fig.l  which  shows  that  the  height  variation  of 
this  parameter  is  double  valued  about  a maximum  at  around  200  km, 
a shape  which  has  important  consequences  for  the  accuracy  of  the 
Faraday  and  protonospheric  content  measurements  which  will  be 
discussed  in  detail  later. 

The  parameter  is  closely  related  to  the  so-called  M-factor 
used  by  many  workers  in  the  reduction  of  Faraday  rotation  angles 


to  equivalent  vertical  electron  content.  A so-called  mean 
ionospheric  height  of  375  km,  at  which  to  evaluate  the  mean  M-factor, 
was  recommended  by  KERSLEY  and  TAYLOR  (1974)  for  Faraday  rotation 


0-4  0-5  0-6  07  0-8  0-9 


f[_  MHz 


Fig.l.  Height  variatior.  ioncitudinal  f^.-rof re'iuer.cy  for 
ATS-6  to  Abcryst  th  path  to  a heifjht  of  2000  laa. 


most  appropriate  for  work  with  geostationary  satellites  altiiougti 
ALMiaOA  (1973)  suggested  350  km. 

Several  workers  (e.g.  MENDILLO  et  al.,  1974  and  P.’J'AGIANNIS 
et  al.,  1975)  have  considered  the  effects  of  a diurnal  variation 
in  the  mean  ionospheric  heiglit  on  measurements  of  Faraday  content. 
This  point  will  be  considered  in  detail  for  tlie  ATS-6  measurements 
later. 

A fixed  value  of  = 0.817  MHz  corresponding  to  a mean 

ionospheric  height  of  420  km  has  been  used  in  the  present  work. 

Thus,  from  Equation  6,  the  relationship  between  the  measured 

polarisation  rotation  (in  degrees)  of  the  140.056  MHz  carrier  and 

-2 

the  Faraday  electron  content  (in  m ) becomes 

N„  = 4.96  X lO'*  n 7. 

F 

(ii)  Modulation  phase  and  total  electron  content 

The  other  angular  parameter  measured  by  the  Aberystwyth 
receiver  is  the  phase  of  a modulation  on  a VHF  carrier  with 
respect  of  that  of  an  identical  modulation  on  a coherent  UHF 
carrier.  After  calibration  to  eliminate  non-ionospheric 
contributions  to  the  measured  phase,  the  ionospheric  phase  delay 
is  given  by 


(j)  _ (^2  - (jlj  - ((Jij^  — (jig)  8. 

where  the  corresponding  suffixes  refer  to  the  frequencies 
f,  = 140.056  iniz,  fa  = 141.0564  MHz,  fa  = 360.144  MHz  and 
f*  = 361.1444  MHz. 

From  Equation  5 the  modulation  phase  delay  for  ordinary 
wave  components  measured  at  Aberystwyth  becomes 
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] Nds 
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. . V r !■  ^ L~  1 + 1 

‘P  - + fi+fL 

The  longitudinal  gyrof requency  is  always  less  than  0.83  Mllz 
for  the  ATS-6  to  Aberystwyth  geometry  and  is  thus  small  compared 
to  the  signal  frequencies.  Hence  if  an  average  value  f is 

la 

chosen  for  the  path  the  above  equation  can  be  written 

(j>  = - k C — — - — ^7-  - — ^ ^ ] / Nds  10. 

fj+fj^  fx+fj^  f,+fj^  fa+fj^  ® 

that  is 

(()  = - k [ — ~ - ^ ^ Nt 

fa+f^  fi+f,  ffc+f^  fa+f, 

u L L 

where  N^  is  the  total  electron  content  along  the  slant  path  from 

the  satellite  to  the  receiving  antennas. 

Using  f.  = 0.817  MHz,  corresponding  to  a mean  ionospheric 

Xj 

^2 

height  of  420  km,  the  total  electron  content  (in  m ) can  be 
obtained  from  the  modulation  phase  delay  (in  degrees) 

N.J,  = - 4.87  X 10'“  4>  12. 

The  mean  value  of  f^^  up  to  2000  km  has  been  estimated  to  be 
0.66  MHz,  however  the  resulting  change  in  the  numerical  constant 
in  Equation  12  is  less  than  0.25%  when  this  value  is  used. 

(iii)  F-f actor 

An  additional  experimental  parameter  which  can  be  obtained 
from  the  ATS-6  angular  data  is  the  shape  factor  or  F-factor  which 
is  a weighted  average  electron  gyrofrequency  along  the  path.  This 
has  been  defined  by  DAVIES  et  al.  (1975a)  as 

= / N f,  ds  / / Nds 

S 1j  s 


F 


13. 


and  is  closely  related  to  the  M-  factor  used  by  SMITH  (1970)  and 
KERSLEY  and  SAMBRUOK  (1971)  which  represents  conditions  along  the 
equivalent  vertical  rather  than  the  slant  path. 

By  substitution  from  Equations  3 and  10  F can  be  related 
to  the  measured  parameters  by 

F = 0.832  5 14. 

9 

Protonospheric  content 

It  can  be  seen  that  absolute  measurements  of  the  total  and 
Faraday  electron  contents  will,  in  principle,  allow  the  electron 
contribution  from  'v  2000  kni  to  the  geostationary  satellite  to  be 
determined.  The  effective  upper  boundary  for  this  region,  the 
protonosphere  or  exosphere,  is  the  plasmapause  which  is  usually 
well  below  the  geostationary  orbit  height  of  L'v6.6.  The  protono- 
spheric electron  content  (Np)  can  be  found  from 

Np  = - Np  15. 

4.  RECEIVING  EQUIPMENT 

a.  Antennas 

The  antenna  system  constructed  at  Aberystwyth  consisted  of 
two  single  element  short-backfire  antennas  for  the  UHF  and  VHF 
transmissions  from  ATS-6.  The  short-backfire  antenna,  originally 
designed  for  microwave  frequencies  (EHRENSPECK  and  STROM,  1971) 

was  chosen  for  the  ATS-6  experiment  because  of  its  wide  bandwidth,  ] 

high  gain,  low  sidelobe  levels  and,  very  importantly  for  modulation  j 

phase  comparisons,  small  group  delay  (GRUBB,  1972).  The  actual  j 

antennas  constructed  followed  the  optimised  designs  of  HARTMANN  5 

and  ENGELHARDT  (1974)  for  the  140  MHx  and  360  MHz  frequencies.  1 


The  circular  framework  for  the  reflectors  was  of  steel  conduit 
covered  with  galvanised  steel  wire  weldmesh  with  a square  mesh 
size  of  2.5  cm.  The.  small  forward  reflector  was  of  sheet  aluminium 
in  each  case  and  the  crossed  dipoles,  also  aluminium,  followed  the 
dimensions  specified  by  Hartmann  and  Engelhardt  and  were  positioned 
to  give  an  antenna  impedance  of  50fJ.  The  outputs  from  the  crossed 
dipoles  were  fed  directly  to  50f2  miniature  balanced  to  unbalanced 
transformers  and  then  combined  in  a quadrative  hybrid  mounted 
inside  the  support  tube  for  the  dipoles  and  reflector.  From  the 
hybrid  50Q  coaxial  feeds  for  both  circular  polarisations  were  taken 
to  the  preamplifiers  mounted  in  a weatherproof  steel  box  at  the 
base  of  the  antennas,  only  the  left  hand  circular  ordinary  mode 
being  used  at  360  MHz.  A photograph  of  the  actual  antennas,  mounted 
on  a small  slope  in  the  satellite  azimuth  at  a viewing  elevation 
of  20°  is  shown  in  Plate  1. 

Approximate  estimates  of  the  antenna  beamwidths  can  be  made  from 
Fig. 2 which  shows  relative  changes  in  received  signal  levels 
averaged  over  24  hours  as  the  satellite  was  moved  through  the 
antenna  beams.  The  satellite  was  started  on  a westwards  drift  on 
2 August  1974  with  a second  manoeuvre  on  4 August  1974  with 
consequent  changes  in  viewing  azimuth  and  elevation.  Fig. 2 indicates 
that,  neglecting  the  small  change  in  elevation,  the  half  beamwidth 
for  a 3dB  drop  in  signal  level  corresponds  to  an  azimuth  change  of 
less  than  15°  for  both  antennas,  a result  which  is  in  broad 
agreement  with  the  estimates  quoted  by  HARTMANN  and  ENGELHARDT  (1974) 
for  short-backfire  antennas. 

This  antenna  beamwidth  suggests  that  for  a viewing  elevation 
of  20°  a ground  reflected  component  of  appreciable  magnitude  may 
be  present  to  interfere  with  the  direct  wave.  The  effects  of  ground 
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27  July  to  23  i\v.ga3t  1976.  The  AT3-6  Jatellito  uas  subject  to  positional  drift 
manocuvers  on  2 and  4 July  1976. 


reflected  signals  on  the  ATS-6  experiment  have  been  discussed  in 
detail  Ly  REES  et  al.  (1977). 


b.  Receiver 


The  design  of  the  receiving  equipment  used  at  Aberystwyth  is 
based  on  the  ATS-6  receiving  system  developed  at  the  Space  Environment 
Laboratory,  NOAA,  Boulder,  USA  (GRUBB  et  al.  1972).  The  actual 
receiver  used  was  constructed  in  modular  form  by  Dunford  Hatfield  Ltd., 
Sheffield,  U.K.  under  the  supervision  of  Mr.  L.  J.  Martin,  Appleton 
Laboratory,  Slough,  U.K. 

The  receiver  uses  phase  locked  loop  techniques  for  the  accurate 

determination  of  phase  differences  between  very  low  level  signals; 

. . - 1 6 — ^ 
the  ATS-6  signal  levels  at  the  ground  being  typically  ^0.6x10  Ita 

(DAVIES  et  al.,  1972).  A schematic  block  diagram  of  the  receiver 

is  shown  in  Fig.  3.  The  Faraday  rotation  and  modulation  phase 

systems  can  be  considered  separately. 


(i)  Faraday  rotation 


In  the  Faraday  rotation  receiving  system  the  two  circularly 
polarised  components  of  the  140  MHz  signal  are  converted  to  an 
intermediate  frequency  of  11.700  MHz  using  a voltage  controlled 
local  oscillator.  The  signal  in  one  i.f.  channel  is  locked  to  a 
reference  fixed  frequency  crystal  oscillator  with  the  out  of 
balance  voltage  from  the  phase  detector  fed  back  with  suitable 
time  constant  to  the  voltage  controlled  local  oscillator.  The 
other  signal  with  opposite  polarisation  is  brought  to  the  same 
i.f.  using  this  common  local  oscillator  and  its  phase  compared  to 
that  of  the  11.700  MHz  reference  in  a phase  detector.  The 
amplified  outputs  in  the  form  of  two  quasi-d.c.  voltages  give  the 
quadrature  components  (X  and  Y)  from  which  the  amplitude  of  the 
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second  signal  (A)  nnd  its  phase  with  respect  to  that  of  the  other  component  (0) 
afe  obtained  from  A = (X^  + Y^)‘  and  0 = arctan  (Y/X) ; the  latter 
yielding  t.lie  Faraday  rotation  (.Q  - 6/2). 

( i i ) Modulation  phase 

The  carriers  and  their  respective  upper  sidebands  are  received 
in  the  modulation  phase  circuits.  The  360  MHz  carrier  and  its 


I sideband  are  brought  to  intermediate  frequencies  of  11.700  MHz  and 

10.700  MHz  respectively  by  common  two  stage  conversion,  the  second 
stage  by  mixing  with  a signal  from  a voltage  controlled  crystal 
oscillator  operating  in  a phase  locked  loop.  The  conversion  of 
one  140  MHz  carrier  to  11.700  MHz  described  above  also  brings  the 
modulation  on  this  carrier  to  an  i.f.  of  10.700  MHz.  The  361  MHz 
sideband,  converted  to  10.700  MHz  is  now  used  in  a loop  to  lock  a 
voltage  controlled  crystal  oscillator  at  this  frequency  which  in 
turn  provides  a reference  for  phase  comparison  with  the  signal  in 
the  i.f.  channel  corresponding  to  the  sideband  on  the  lower  frequency 
carrier.  In  this  way  the  quadrature  component  outputs  from  a 
phase  detector  give,  the  phase  of  the  1 MHz  modulation  on  140  MHz 
I with  respect  to  that  on  360  MHz. 

The  entire  receiving  system  was  designed  as  a modular  integrated 
package  in  which  the  circuits  have  been  simplified  by  using  common 
modules  and  avoiding  duplication  of  channels. 

The  original  equipment  allowed  for  acquisition  of  the  signals 
by  manual  opening  of  the  phase  locked  loops  and  adjustment  of 
suitable  d.c.  levels  fed  to  the  summing  amplifiers.  However, 
additional  circuitry  was  built  for  automatic  signal  acquisition 
allowing  the  equipment  to  be  used  reliably  at  an  unattended  field 
station  with  the  often  interrupted  ATS-6  transmissions. 
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Thft  settings  of  tlie  time  constants  in  the  p)iase  locked  loops 
are  a compromise  between  narrow  loop  bandwidth  and  thus  improved 
signal  to  noise  ratio,  and  rapid  signal  acquisition  and  holding. 

It  was  found  tliat  loop  bandwidth  settings  of  3 Hz  for  the  two 
carrier  loops  and  10  Hz  for  tlie  sideband  loop  represented  the 
optimum  compromise  in  the  operation  of  the  receiver.  The  attenuators 
in  the  i.f.  channels  were  adjusted  at  an  early  stage  to  give  signal 
amplitudes  at  the  outputs  in  the  range  3 to  5V  and  thereafter  kept 
constant.  The  quadrature  outputs  from  the  phase  detectors  are 
brought  to  a typical  level  of  several  volts  in  the  d.c.  amplifiers. 

In  the  original  design  the  amplifier  time  constant  was  Is,  however 
this  was  modified  to  10s  for  the  phase  measurements  in  the 
Aberystwyth  system  providing  some  smoothing  of  the  output  voltages. 

c . Data  Logging,  Coordinate  Converters  and  Recorders 

The  output  levels  of  several  volts  are  convenient  for  storage 
of  data  after  analogue  to  digital  conversion  on  magnetic  tape  for 
subsequent  computer  processing.  Readings  of  six  voltages  corres- 
ponding to  two  pairs  of  quadrature  components  and  two  carrier  levels 
were  stored  on  the  data  logging  system.  For  quick  look  monitoring 
of  the  performance  of  the  equipment  and  ionospheric  changes  each 
pair  of  quadrature  components  were  fed  to  an  analogue  rectangular 
to  polar  coordinate  converter  vjhose  output  voltage,  a linear 
function  of  angle,  was  fed  to  a pen  recorder.  A two  channel 
Scrvoscribe.  chart  recorder  was  used  to  give  simultaneous  records 
of  Faraday  rotation  and  modulation  phase  angles.  A chart  width 
of  20  cm  for  both  channels,  corresponding  to  360°  phase  difference. 


gave  good  angular  resolution  while  the  temporal  speed  was  2mm  min 
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Sariplc  chart  rccordinc  of  noclaJati  on  phaoe  and  Faraday  rotation 


A sample  record  is  shown  in  Fig.  A,  A second  recorder  monitored 
tiie  amplitudes  of  the  140  MHz  and  360  M}iz  carriers. 


5.  CALIBRATION  AND  ACCURACY 

In  the  ATS-6  experiment  the  accuracy  of  the  absolute  measurements 
of  modulation  phase  and  Faraday  rotation  is  of  great  importance.  As 
the  protonospheric  electron  content  is  determined  from  the  difference 
between  the  total  and  Faraday  contents,  and  during  daytime  may 
represent  only  10%  to  20%  of  the  former,  to  achieve  an  accuracy 
of  10%  in  N requires  the  determination  of  N’  and  N„  to  about  1%. 

P if 

The  overall  accuracy  of  ATS-6  measurements  is  dependent  on  various 
factors  like  the  spacecraft  transmitter  and  antennas,  the  performance 
of  the  receiving  antennas  and  the  calibration  of  the  receiving 
system.  These  will  be  discussed  below.  In  addition,  the  effect 
of  the  choice  of  mean  ionospheric  height  on  the  Faraday  observation 
will  be  considered  in  detail  in  the  interpretation  of  the  results. 

The  calibration  and  accuracy  of  the  modulation  phase  and  Faraday 
rotation  measurements  are  considered  separately  below. 

a.  Modulation  Phase 

The  ionospheric  phase  delay  (<j>)  of  the  141  MHz  signal  can  be 
obtained  from  the  phase  angle  measured  at  the  receiver  outputs 

if  various  systems  calibrations  are  known.  The  relationship 
can  be  written  in  the  form 

<)>  = - Up  + ‘t>A  'i’c  '*'r^  “ 360  X M 16 


where  ij)^  represents  the  transmitter  delay,  takes  account  of 
the  group  delays  of  the  receiving  antennas  and  their  physical 
separation,  <(i  allov;s  for  differences  in  cable  lengths,  gives 


the  receiver  delay  and  the  final  term  with  M integer  accounts 
for  complete  cycle  ambiguities.  Each  of  these  factors  can  be 
treated  separately. 

(i)  Transmitter  delay 

The  ATS-6  radio  beacon  was  subject  to  careful  prelaunch 
calibration  to  ascertain  the  phase  relationship  of  the  signals 
at  the  transmitting  antennas  and  the  long  term  phase  stability 
of  the  transmitter.  The  transmitter  delay  was  found  to  be 
41^  = - 68°  (DAVIES  et  al.,  1975b)  with  a temperature  stability 
of  the  ± 2°  over  the  expected  operating  range.  This  value  of 
transmitter  delay  has  been  used  in  the  present  study. 

(ii)  Receiving  antennas 

The  choice  of  short-backfire  antennas  for  the  receiving 
system  was  conditioned  in  part  by  the  favourable  phase/frequency 
characteristic  of  this  design  in  comparison  to  other  high  gain 
antenna  systems. 

The  problem  of  antenna  phase  centre  displacement  for  carrier 
and  sideband,  representing  a frequency  dependence'  of  s in  the 
first  term  of  Equation  A,  and  its  interpretation  in  terms  of 
antenna  group  delay  has  been  discussed  in  detail  by  GRUBB  (1975). 
Experimental  measurements  of  the  group  delays  of  short-backfire 
antennas  of  optimised  design  for  the  ATS-6  frequencies  have  been 
made  by  HARTMANN  and  ENGELHARDT  (1974  and  1975).  They  report 
a group  delay  of  5.8  ns  for  the  1 MHz  sideband  on  360  MHz  while 
the  corresponding  delay  for  the  140  MHz  antenna  was  15  ns,  values 
which  scale  approximately  according  to  wavelength.  Thus  the 
grou])  delay  appropriate  to  ATS-6  measurements  has  a magnitude 
15  - 5.8  = 9.2  ns,  which  for  a modulation  frequency  of  1 MHz 
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amounts  to  a phase  delay  of  -2.9°. 

The  antennas  for  360  Mllz  and  140  MHz  were  mounted  with  their 
active  elements  coplanar  normal  to  the  propagation  direction. 

Since  the  phase  centre  of  the  short-backfire  can  be  considered  to 
be  in  the  plane  of  the  large,  reflector  (GRUBB,  1975)  the  Aberystwyth 
antennas  must  be  considered  as  physically  separate  in  the  propagation 
direction  with  the  360  MHz  antenna  ahead  of  the  other  by  'v  0.34m. 

This  corresponds  to  an  additional  delay  1 ns  which  represents 
a phase  at  1 MHz  'v  - 0.3°.  Thus  to  an  accuracy  appropriate  to 
the  experiment,  allowing  for  neglect  of  the  second  order  term 
evaluated  for  only  one  frequency  by  HARTMANN  and  ENGELllARDT  (1974) 
and  taking  into  account  small  differences  in  antenna  design  and 
feed  arrangements  compared  to  those  of  Hartmann  and  Engelhardt, 
the  value  of  used  in  the  present  study  has  been  taken  to  be  -3°. 

As  identical  lengths  of  cables  were  used  in  comparable 
situations  throughout  the.  receiving  system  (jj,  has  been  taken  to 
be  zero  in  the  interpretation  of  the  Aberystxvyth  data. 

(iii)  Receiver  calibration 

To  obtain  the  group  delay  of  the  receiving  system,  from  the 
inputs  to  the  preamplifiers  to  the  outputs  to  the  data  logger  and 
chart  recorder,  it  is  necessary  to  generate  a spectrum  of  signals 
to  simulate  the  satellite  transmissions.  This  is  done  by  a 
receiver  calibrator  which  was  built  as  part  of  the  receiving 
equipment  to  the  specifications  of  the  design  of  GRUBB  €>.  al.  (1972). 

In  use,  the  calibrator,  kept  physically  well  away  from  the 
actual  receiver  to  minimise  cross  coupling  of  stray  signals,  was 
connected  to  the  inputs  to  the  preamplifiers  via  the  attenuator 
and  power  splitter  network  shown  in  Fig.  5a.  The  attenuators 
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were  adjusted  to  give  receiver  output  signal  levels  for  the  two 
carriers  comparable  to  those  obtained  with  the  actual  satellite 
transmission  and  the  resulting  modulation  pliase  delay  for  the 
receiver  determined.  Calibrations  carried  out  througliout  the 
entire  period  of  ATS-6  observation  gave  consistent  results. 
Additional  comparisons  were  made  on  several  occasions  using  an 
ALDI  design  calibrator  on  loan  from  Air  Force  Geophysical 
Laboratory,  Bedford,  USA  which  in  general  gave  consistent  results 
with  the  Aberystwyth  calibrator  and  discrepancies  could  be  traced 
to  the  performance  of  the  ALDI  instrument.  Tests  with  the  master 
standard  calibrator,  made  available  to  ATS-6  observers  by  the 
Boulder  group,  gave  a modulation  phase  delay  for  the  receiver 
which  was  within  the  error  limits  obtained  using  the  Aberystwyth 
calibrator. 

It  was  noticed  in  the  operation  of  the  ATS-6  receiver  that 
amplitude  and  phase  were  not  strictly  independent  in  that  the 
modulation  phase  angle  showed  a small  dependence  on  the  360  MHz 
carrier  amplitude.  The  form  of  this  dependence,  which  appeared 
to  be  attributable  to  the  360  MHz  preamplifier,  can  be  seen  from 
Fig. 6 which  shows  the  receiver  calibration  modulation  phase  angle 
as  a function  of  signal  level  defined  by  the  setting  of  the 
attenuator  in  the  360  MHz  arm  (Fig.  5a).  The  results  displayed 
in  Fig. 6 were  obtained  using  the  standard  Boulder  calibrator. 

The  individual  points  give  a mean  value  and  its  standard  deviation 
obtained  from  the  average  of  some  60  readings  taken  over  a period 
of  about  5 minutes.  Comparable  data  for  the  variation  of  signal 
level  in  the  lAO  MHz  channel  presented  in  Fig. 7 shows  an  essential 
independence  of  phase  and  amplitude. 


360  MMz  Attenuation  (t!B) 


Pif;.6.  Roccivor  modulation  phaco  delay  as  a fua.otion 
of  attenuation  in  360]'I!b  channel. 
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140  MHz  Attenuation  IdB)  !j 


Fif;.7.  Roceivor  no-^ulation  Tihaac  delay  as  a function  of 
attenuation  in  KORIIe  channel,  with  0 ~ 10  (IB 
attenuator  (Fit';. 5a)  in  I4OF1H2  am, 
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The  received  signal  level  at  360  MHz  was  subject  to 
variations  within  a typical  range  of  about  3dB  which  have  been 
attributed  to  changes  in  satellite  pointing  and  possibly  radiated 
power,  so  that  a small  error  in  modulation  phase  is  introduced 
from  the  amplitude/phase  coupling,  labile  it  is  possible  to  take, 
account  of  this  variation,  for  the  purposes  of  the  present  study 
which  is  primarily  concerned  with  averaged  observations,  a 
calibration  value  has  been  used  which  is  representative  of  the 
approximate  mean  360  MHz  signal  level  received  throughout  tlie 
observing  period. 

In  the  long  term,  the  receiver  calibration  was  found  to  have 
good  stability  with  only  small  dependence  on  external  temperature. 

Taking  account  of  the  factors  discussed  above  a value  of 

modulation  phase  delay  in  the  receiver  <}i„  = ~ has  been  used 
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for  the  analysis  of  the  results  described  in  this  report.  Tt  is 
believed  that  this  value  is  representative  of  the  receiver 
calibration  applicable  to  this  study  of  averaged  data  throughout 
the  nine  months  of  observation  to  an  accuracy  of  about  ± 3°. 
However,  in  the  investigation  of  isolated  individual  events, 
particularly  at  times  of  scintillation  of  the  360  MHz  signal, 
the  errors  could  be  somewhat  larger. 


b.  Faraday  Rotation 


The  calibration  of  the  Faraday  rotation  system  to  obtain 
absolute  values  of  measured  phase  difference  of  the  ordinary  and 
extraordinary  modes  involves  consideration  of  the  orientations 
of  the  transmitting  and  receiving  antennas  together  with  the  phase 
delay  of  the  receiver. 


(i)  Transmitting  antenna  orientation 


The  me-surcment  of  the  radiated  polarisation  was  carried  out 
as  one  of  the  pre-launch  calibrations  of  the  satellite  radio  beacon 
so  that  its  orientation  is  kno\im  with  respect  to  the  satellite 
axes  (DAVIES  ct  al.  1975a).  It  is  thus  possible  to  calculate  the 
orientation  of  the  transmitted  polarisation  as  viewed  by  an  observer 
for  known  satellite  and  observer  locations.  For  geometry  applicable 
to  the  ATS-6  satellite  stationed  at  35^E  viewed  from  Aberystwyth 
the  effective  140  MKz  transmitting  antenna  orientation  makes  an 
angle  of  26°  with  the  local  station  horizontal  measured  in  a counter 
clockwise  direction  looking  towards  the  satellite. 

( i i ) Receiving  antenna  orientation 

The  active  elements  of  the  140  MHz  receiving  antenna  are  crossed 
linear  dipoles  followed  by  a quadrature  hybrid  in  which  the  linearly 
polarised  incoming  signal  is  converted  to  its  circularly  polarised 
components.  A zero  phase  difference  between  these  polarisation 
modes  implies  that  the  angle  of  the  received  linear  polarisation 
bisects  the  crossed  dipoles.  The  dipoles  were  set  to  be  symmetrically 
inclined  at  nominal  angles  of  ± 45°  with  respect  to  the  local 
horizontal,  however  measurements  with  an  inclinometer  showed  that 
there  was  an  assymmetry  of  2°  in  the  adjustment.  Thus  the  angle 
corresponding  to  zero  phase  difference  between  the  circular  modes 
was  found  to  be  92°  from  the  local  horizontal,  counter-clockwise 
viev;ed  from  behind  the  antenna,  that  is  in  the  direction  of 
increasing  Faraday  rotation. 

Thus  the  difference  between  the  effective  polarisations  of 
the  transmitting  and  receiving  antennas  amounts  to  92  - 26  = 66° 
which  must  be  allov;cd  for  in  estimating  absolute  Faraday  rotation. 
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(iii)  Receiver  calibration 


The  phase  comparison  of  the  signals  in  the  140  Mllz  channels 
of  the  receiver  was  calibrated  using  the  calibrator  described 
earlier.  A signal  from  the  calibrator,  brought  to  an  appropriate 
level  in  attenuators  was  fed  from  a power  splitter  to  the  inputs 
of  the  two  140  MHz  preamplifiers  (Fig.  5b)  and  the  corresponding 
output  from  the  Faraday  rotation  channels  noted.  A value  of  342° 
phase  difference  between  the  polarisation  modes  corresponding  to 
a Faraday  rotation  angle  of  171°  has  been  used  for  the  present 
study.  This  value  was  representative  of  the  early  calibration 
runs  to  an  accuracy  of  about  ± 2°  rotation.  However,  it  should 
be  noted  that  angles  outside  this  range  towards  lov;er  values  were 
obtained  in  later  calibration  work  with  an  extreme  value  of  166° 
being  measured  on  one  occasion.  The  physical  situation  of  the 
receiving  equipment  at  Aberyst\jyth  in  a small  field  site  hut  made 
regular  calibration  impractical  so  that  it  has  been  necessary  to 
assume  a fixed  value  for  the  present  work. 

An  additional  factor  affecting  the  polarisation  measurement 
was  the  operation  for  some  months  of  a second  group  delay  receiver 
of  ALDI  design  in  parallel  with  our  own  equipment  using  the 
extraordinary  mode  polarisation  from  the  same  antenna  system.  It 
was  found  that  a change  of  about  2°  was  introduced  in  the 
polarisation  angle  measurement  when  one  of  the  140  MHz  signals 
was  divided  between  the  two  receivers. 

An  independent  check  of  the  Faraday  rotation  calibration  was 
made  using  a separate  receiving  system  fed  from  a dipole  which 
could  be  manually  rotated.  The  polarisation  angle  of  the  radiation 
received  from  the  satellite  was  determined  by  swinging  the  dipole 
through  a signal  null  between  equal  levels  at  the  receiver  output. 
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Measurements  over  an  extended  period  with  changing  polarisation 
were  then  compared  with  the  chart  record  of  Varaday  rotation. 

The  results  of  one  such  calibration  check  representing  some  80 
independent  measurements  taken  over  a 10  hour  period  are  plotted 
in  Fig.  8.  A best  fit  straight  line  yields  an  intercept  of 
13°  ± 2°  for  the  angular  difference  between  the  two  measurements. 
This  angle  can  be  compared  to  a corresponding  value  of  11°  obtained 
from  the  calibration  measurements  of  receiver  and  receiving  antenna. 

c.  Resolution  of  Ambiguities  • 

Having  determined  the  origins  for  the  measurements  of 
modulation  phase  and  polarisation  rotation  one  further  consideration 
remains  in  the  determination  of  absolute  values  of  these  angles, 
namely  ambiguities  of  complete  cycles.  This  point  dictated  the 
choice  of  frequencies  for  ATS-6  and  its  resolution  was  facilitated 
by  the  recording  of  two  angles  iv^ith  differing  sensitivities  to 
ionospheric  changes  so  that  it  was  relatively  easy  to  select  the 
correct  cycles  to  give  reasonable  values  for  both  ionospheric  and 
protonospheric  parameters.  Furthermore,  the  results  reported  here 
were  obtained  for  conditions  of  minimum  solar  activity  so  that 
the  predawn  modulation  phase  angles  were  usually  within  the  first 
cycle.  Thus  in  spite  of  the  frequent  interruptions  of  the 
transmissions  the  resolution  of  complete  cycle  ambiguities  has 
presented  no  real  problems. 

6 . IONOSPHERIC  PARAMFTF.RS  - RESULTS 

a.  Data  Base 

The  ionospheric  data  presented  here  have  been  obtained  from 
the  chart  records  of  the  angular  parameters  ()i  and  0 reduced  at  one 
hour  intervals  for  the  period  November  1975  to  July  1976  v;hen 


Receiver  A.' 


Fi^.8.  AT3-6  receiver  Faraday  rotation  nnfrles  compared 

to  independent  raoasurenents  obtained  unin.'c  a 
rotating  dipolo. 


satellite  transmissions  were  available  and  records  were  obtained. 


Hourly  values  of  N^,  Nj^,  and  together  with  F have  been 

calculated.  The.  F-f actor  has  been  presented  in  two  equivalent 

forms,  F itself  and  h_,  (or  h rf  after  DAVIES  et  al . (1976)),  the 

dL  eti 

height  for  which  the  observed  value  of  F equals  the  longitudinal 
gyrof requency . Tables  of  these  parameters  on  a monthly  basis 
have  been  prepared  which  also  contain  a statistical  summary  of 
the  data  for  each  hour,  namely  mean  and  standard  deviation, 
median  and  quartile  values.  Sample  tabulations  of  N , N_,  N and 
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F for  the  month  of  November  1975  are  presented  in  Tables  2 to  5. 
b . Day  to  Day  Variability 


The  day  to  day  variability  of  the  parameters  can  be  gauged 

from  superposed  multiplots  of  the  diurnal  variation  for  each  day 

in  a month.  Plots  of  N„,  N'  and  N for  three  months,  representing 
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winter,  equinox  and  summer  are  shown  in  Fig.  9.  The  large  day  to 
day  variability  of  N and  N„  agrees  with  the  findings  of  many 
beacon  satellite  workers  but  as  yet  no  completely  satisfactory 
explanation  is  available  (RISHBETH  and  KOHL,  1976).  Apart  from 
obvious  solar  and  geomagnetic  dependencies  it  has  been  suggested 
that  the  variability  reflects  the  erratic  nature  of  the  global 
neutral  wind  system  (KANE,  1975). 

Fig.  9 shows  that  the  variability  is  most  pronounced  in  Np, 
to  an  extent  that  the  general  shape  of  the  individual  diurnal 
plots  is  masked  by  the  day  to  day  variations.  In  addition,  the 
local  mean  times  of  the  maximum  and  minimum  values  of  Np  show 
considerable  ranges,  although  a more  consistent  picture  is  found 
in  the  monthly  averaged  data. 
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ilultiplots  o-T  diumnl  vari:itionc  of  N. 
Abory.'-t’.pv'th  for  noloctod  noni.hs. 


In  Table  6 standard  deviation  of  tVic  noon  values  of  each  of 


the  parameters  by  month  is  compared  to  the  corresponding  diurnal 
range  of  the  monthly  mean  values.  For  N and  N the  variability 
is  much  less  than  the  diurnal  range  with  the  smallest  standard 
deviations  in  winter  and  largest  at  equinox  suggesting  that  the 
magnitude  of  the  day  to  day  variations  is  related  to  the  electron 
content  level.  However,  for  Np  in  winter  the  noon  time 
variability  is  comparable  to  the  mean  diurnal  range  with  the 
deviations  representing  some  30Z  to  50%  of  this  range  at  other 
seasons.  The  variability  of  Np  from  day  to  day  is  dependent  on 
many  factors  and  in  discussing  this  parameter  account  must  be 
taken  of  the  local  and  conjugate  ionospheres  together  with  the 
rapid  depletion  and  slow  refilling  of  the  plasmasphere  after 
geomagnetic  activity,  points  which  will  be  considered  in  detail 
later.  In  addition,  double  peaked  diurnal  curvcsand  nighttime 
maxima  in  Np  are  observed  on  occasion,  features  which  correspond 
to  the  light-ion  spectrometer  measurements  of  CHEN  et  al . (1976), 
however  such  variations  are  smoothed  out  in  the  median  curves. 

c . Monthly  Median  Behaviour 

Because  of  the  large  day  to  day  variability  in  the  parameters 
it  was  decided  to  concentrate  study,  at  least  initially,  on 
values  representing  averaged  data  for  each  month. 

The  diurnal  plots  of  the  median  values  for  each  month  are 

plotted  in  Figs.  10  to  18.  These  show  N N N , N /N  %,  F and 

1 , i*  , i r 1 

h_-  (i.e.  h for  each  of  the  nine  months  under  study.  The 
BL  eff 

dotted  sections  of  the  curves,  particularly  for  March  1976,  join 
median  points  based  on  less  than  ten  individual  daily  values. 


Table  6 


Standard  deviations  of  the  noon  values  of 


N.p,  Nj,  and  N^,  compared  to  diurnal  range  of 


mean 

data  for 

each 

month. 

N,j,  10 

1 

-m  1 

i 

N„  10‘*  m"* 

Np  10 

RMS 

Range  ' 

RMS 

Range 

RMS 

Range 

Nov.  75 

3.3 

16.5  1 

2.8 

15.8 

7.5 

8.3 

Dec.  75 

2.4 
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10.7  1 

1.7 

10.3 
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9.5 

j Jan.  76 

2.0 

1 
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10.0 
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2.2 
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2.2 
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5.2 
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3.1 
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3.2 
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5.3 

14.8  1 

j Apr.  76 

4.1 
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3.5 

19.0 

9.5 
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1 May  76 
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while  sections  missing  iiulicate  that,  no  transmissions  were  availal)le 
at  these  hours  throughout  the  month. 

The  diurnal  variations  of  and  are  essentially  similar 
in  character  indicating  that  as  expected  the.  ionospheric  contribution 
dominates  the  total  electron  content  at  least  during  daytime.  By 
contrast,  Np  shovv’s  only  a small  diurnal  variation  with  minimuni 
values  in  the  predawn  hours  around  3 x 10  m for  all  months 

— 2 

and  a maximum  generally  in  late  afternoon,  of  about  A x 10 m 

_ 2 

in  winter  rising  to  approximately  5 x 10^®  m in  summer.  The 
protonospheric  content  expressed  as  a percentage  of  the  total 
has  values  by  day  typically  in  the  range  15  to  20%  while  the 
predavm  maximum  of  about  35%  in  summer  exceeds  40%  in  winter. 

The  diurnal  variations  of  F show  a predawn  minimum  about 
0.5MHz  corresponding  to  an  effective  height  of  about  1800  km  while 
the  daytime  value  is  generally  in  the  region  of  0.66  Milz  repre- 
senting an  effective  lieight  around  1200  km.  It  is  interesting 
to  note  that  F is  approximately  proportional  to  the  factor 
studied  by  KERSLEY  and  SAMBROOK  (1971)  who  found  a 30%  diurnal 
variation  in  this  parameter  wliich  is  in  close  agreement  with  the 
present  measured  change  in  F.  The  effective  heights  can  be 
compared  to  those  quoted  by  HARGREAVES  and  HOLMAN  (1972)  from 
model  studies  which  were  in  the  range  1000  to  1400  kin. 

d.  Seasonal  Variations 

The  daytime  ''alucs  of  N„  and  N„  are  smallest  in  winter  with 
the  maximum  just  after  noon.  The  larger  single  peak  later  in 
the  afternoon  at  equinox  gives  way  to  a twin  peak  diurnal  variation 
in  summer  when  the  polewards  neutral  wind  causes  a depression  in 
the  early  afternoon  values.  The  seasonal  variations  of  ionospheric. 

N,,  are  best  appreciated  by  tlie  contour  plot  of  the  median  values 
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ils  a liiiirtioii  ol'  time  of  day  and  nioiilh.  Superposed  on  Fig.  19 
are  the  limes  of  sunrise  and  sunset  at  the  300  km  sub-ionospheric 
point  estimated  from  COLIN  ct  al . (1966).  The  basic  solar 
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control  of  the  N^,  photo-ionisation  is  immediately  apjiarent. 

No  clear  relationship  is  seen  between  each  of  N„,,  K„  and 

i r 

Np  and.  either  mean  solar  flux  or  mean  although  changes  in 
solar  activity  are  small  over  this  nine  month  period  at  solar 
minimum. 

The  noon  time  values  of  shovj  a variation  by  month  which 

follows  closely  the  corresponding  changes  in  N . However  it  is 

r 

apparent  from  an  examination  of  a contour  plot  of  Np  by  hour  and 
month  (Fig.  20)  that  this  parameter  does  not  show  the  same  direct 


solar  control  as  N . Furthermore,  the  daytime  maximum  in  N.,  in 

r i 

Novemiior  and  December  extends  considerably  beyond  the  sunlit 
hours  of  the  short  winter  day,  so  that  even  if  the  time  constants 
involved  are  several  hours  it  is  still  difficult  to  explain  the 
observed  diurnal  variation  purely  in  terms  of  interaction  with 
the  local  ionosphere. 

Several  studies  have  been  concerned  with  plasma  exchange 
between  two  conjugate  ionospheres  along  a field  lino  (c.g.  YOUNG 
et  al.,  1970,  PARK,  1970,  MAYR  ct  al.,  1972)  and  it  is  reasonable 
to  assume  that  the  intervening  protonosphere  is  dependent  on  both 
local  and  conjugate  ionospheric  conditions.  In  a first  attempt 
to  account,  at  least  qualitatively,  for  the  general  form  of  the 
Np  variations,  particularly  for  winter,  ionospheric  data  for 
Slough  (31.5°N,  0.6*'’w)>  representing  conditions  in  the  general 
vicinity  of  the  local  300  km  sub-ionospheric  point  (47.9°N,  2.5”k) 
was  considered  in  conjunction  with  similar  data  for  Cape  Town 


time, and  taken  Co  be  indicative  ot  conditions  in  the  conjugate 
ionosphere  (36°S,  2o'’h) . As  no  current  ionosplieric  data  were 
available  for  Cape  Town  corresponding  values  from  the  previous 
solar  cycle  were  used.  Various  studies  were  made  using  these, 
data  sets  and  it  was  found  that  a contour  plot  which  broadly 
shows  some  of  the  features  of  the  variations  is  obtained 
using  averaged  values  of  corresponding  maximum  electron  densities 
for  the  two  stations  which  are  time  delayed  by  about  two  hours. 

In  particular  the  length  of  the  summer  day  in  the  conjugate 
hemisphere  helps  account  for  the  Np  variations  in  November  and 
December.  The  results  of  this  work  suggest  that  account  must  be 
taken  of  both  local  and  conjugate  ionospheres  in  an  explanation 
of  the  observed  Np  variations.  However,  tliis  study  has  not  been 
pursued  further  in  this  form  because,  concentrating  as  it  does 
on  only  one  field  lino, it  has  limited  application  to  the  actual 
physical  situation  and  furtherm>5rc  its  application  requires  a 
more  realistic  assessr.;ent  of  the  time  constants  involved  in 
protonospher ic  filling  and  depletion. 

e.  Accuracy  of  Np  Data 

Since  N is  obtained  from  the  difference  of  N„  and  N„  and, 

F Jr 

during  daytime  at  least,  has  a magnitude  generally  less  than  20% 
of  N.J,,  it  is  important  to  consider  the  accuracy  of  the  measurements 
of  N„  and  N„  and  the  consequent  accuracy  in  N’  . The  errors  fall 
into  two  categories,  systematic  observational  errors  in  the 
absolute  determination  of  the  angular  parameters  and  errors 
associated  with  the  choice  of  mean  ionospheric  height  in  the 
reduction  of  the  Faraday  rotation  angles. 
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(i)  Observational  errors 
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The  eifects  of  errors  in  the  absolute  values  of  the  measured 
angular  parameters  have  been  estimated  and  their  implications  on 
the  monthly  median  diurnal  variations  of  Np  considered. 

Fig.  21  shows  the  effects  of  various  constant  systematic 
errors  in  the  measured  angles  on  the  diurnal  variations  of  the 
median  Np  values  for  December  1975,  April  1976  and  July  1976. 

Each  curve  represents  the  percentage  error  in  Np  resulting  from 
an  assumed  error  in  one  of  the  angles  iji  and  J2;  the  curves  for 
an  error  of  given  sign  in  (p  and  an  equal  error  of  opposite  sign 
in  52  are  virtually  identical  since  the  numerical  constants  in 
Equations  7 and  12  are  almost  equal.  It  is  believed  that  the 
errors  in  the  measured  parameters  of  the  present  study  can  be 
conservatively  estimated,  in  most  circumstances,  at  less  than 
± 7°  in  each  of  the  angles  so  that  the  resulting  absolute 
accuracy  i>-f  the  Np  data  is  generally  better  than  about  ± I5Z. 

The  short  terra  relative  errors  are,  in  general,  much  smaller  than 
this  value. 

(ii)  Effects  of  constant  f,  factor 

■ ■ ' ' ' ' “ ' ■ ■ 1 ■ - - --  ■ _ ^ -r-.  -T  — 

An  additional  error  could  arise  in  the  data  from  the  use  of 
a constant  f^^  in  the  Faraday  content  estimated.  Two  factors 
must  be  considered  here.  Firstly  there  are  the  consequences  of 
a change  in  the  mean  ionospheric  height  in  the  choice  of  a 
constant  f , and  secondly  the  implications  of  a changing-height 
of  the  F2-layer  peak  with  a resulting  diurnal  variation  in  f 
must  be  'iiscussed.  This  latter  point  is  of  particular  importance 
in  ascertaining  the  validity  of  the  small  diurnal  variations 
found  in  Np. 


Examination  of  Fig.  1 shows  that  in  the  region  300  to  500  km, 
the  value  of  differs  from  the  value  at  420  km  by  less  than 
± 1.3%,  so  that  resulting  uncertainty  in  Np  during  daytime  is  less 
than  10%  and  much  smaller  at  night.  Most  v/orkers  have  chosen  a 
mean  ionospheric  height  between  350  to  400  km,  and  in  this  range 
the  dependence  of  f,  on  height  is  small  for  the  ATS-6  to  Aberystt^yth 
geometry,  representing  a change  in  f^  of  less  than  1%.  The  effect 
on  the  diurnal  variation  of  Np  and  Np/Np%  of  the  use  of  different 
mean  ionospheric  heights  over  a wide  range  from  200  to  600  km  is 
shown  in  Figs.  22  to  24  for  winter  (December,  1975),  equinox 
(April  1976)  and  summer  (July  1976). 

The  height  of  the  maximum  ionisation  in  the  F-layer  shows 
a diurnal  variation  so  that  it  is  reasonable  to  assume  a similar 
change  in  the  mean  ionospheric  height  for  the  precise  evaluation 
of  f in  Faraday  content  measurements.  KERSbEY  and  TAYLOR  (1974) 

Li 

compared  incoherent  scatter  profiles  with  Faraday  data  for 
satellites  in  lovj  orbit  and  found  that  on  average  the  mean  iono- 
spheric height  was  some  80  km  above  the  height  of  the  layer  peak  (h^) 
although  the  difference  between  the  heights  for  individual 
measurements  ranged  from  60  to  120  km.  In  the  present  study  h^ 
has  been  determined  from  ionosonde  data  and  the  consequences  on 
Np  of  the  use  of  three  height  increments,  50,  100  and  150  km 
above  the  layer  peak  assessed. 

To  obtain  monthly  median  representative  values  of  h^,  a 
true  height  reduction  of  the  ionosonde  data  from  Slough  has  been 
performed,  employing  a model  suggested  in  its  original  form  by 
SHIMAZAKI  (1955) . The  virtual  height  of  reflection  is  given  by 


176  km 


t 


h'F2  = 


M(:iOOC))l'2 


where  M(3000)l'’2  is  the  maximum  usable  frequency  factor  in  MHz. 

Tlie  virtual  height  at  a frequency  of  0.834  f^F2  is  numerically 

equal  to  h in  the  absence  of  E and  FI  ionisation.  BR/dlLEY  and 
m 

DUDENEY  (1973)  introduced  a correction  to  account  for  the  daytime 
E-layer  ionisation  and  it  is  this  equation  which  has  been  used 
in  the  present  work.  Their  expression  gives  h^^  directly,  with 
errors  less  than  20  km,  and  is  of  the  form 


. _ 1490  , 

^m  M(3000)F2+AM 


with  AM  = 0.18/(X-1.4)  and  X = f F2/f  E. 

Median  hourly  values  of  h^  on  a monthly  basis  have  been 
obtained  using  Slough  ionosonde  data  and  the  diurnal  variation 
of  f assessed  for  the  three  constant  height  increments.  The 

I4 

resulting  curves  showing  the  effect  of  layer  height  changes  on 
the  diurnal  variation  of  Np  are  presented  in  Figs.  25  to  27  for 
the  three  months  representing  winter,  equinox  and  summer.  These 
show  that  the  diurnal  variation  found  in  Np  is  enhanced  by  the 
use  of  a variable  mean  ionospheric  height  when  compared  to  the 
values  obtained  for  a constant  height.  Following  the  work  of 


KERSLEY  and  TAYLOR  (1974)  an  increment  of  100  km  would  seem 
appropriate  w'hen  using  data  from  a geostationary  satellite  so 
that  the  lowering  of  the  F2-layer  peak  during  daytime  appears 
to  give  values  of  Np  increased  by  some  6%  while  at  night  the 


values  depart  little  from  those  obtained  with  a fixed  f^^ 


evaluated  at  420  km. 

The  important  conclusions  for  this  work  are  that  for  the 
Aberystwyth  observations  the  values  of  Np  are  not  particularly 


'‘P/  ’ ''‘P"‘20 


/••/^  - / 


!,,,,+ 150 

Iiin  -Moo 

h,,!-*-  iO 

420 


Fig. 27.  Kffoct  of  variable  evaluated  for  three  height 
ir.crocioata  to  h , on  riiui'nal  variation  of  I!p  for 
Aberyr.twy  th  for'^July  197&. 


sensitive  to  the  choice  of  mean  i onar.v)heri e height  and  the 
diurnal  variation  oi  this  parameter  is  enhanced  by  the  use  of 
a variable  'leiglit  so  that  it  can  be  concluded  that  this  observed 
diurnal  variation  represents  a true  reflection  of  protonospher ic 
behaviour  and  is  not  simply  a product  of  tlie  observational 
geometry  and  the  analysis  techniques. 


7.  lONOSPilhllTC/rUOTOMOSPHllRfC  FLUX 


a.  Introduction 


Plasma  produced  by  photo-ionisation  in  the  ionospheric  part  of 
the  plasmasphcre  can  travel  up  or  dotm  by  three  possible  mechanisms. 

1.  Electric  fields,  originating  from  either  the  dynamo  region  or 
at  higher  latitudes  the  solar  wind  induced  davat-dusk  field  in  the. 
magnetosphere,  can  cause  vertical  plasma  drifts.  However,  the 
magnitude  is  small  at  mid-l  atitude.s  because  of  the  inclination 

of  the  geomagnetic  field  lines  but  outside  the  plasr.iapause  the 

plasma  can  escape  freely.  EVANS  (1972)  has  reported  electric 

field  induced  vertical  drifts  in  the  F-region  with  a magnitude 
_ 1 

of  about  25  ns 

2.  The  thermospheric  neutral  wind,  resulting  from  atmospheric 
horizontal  temperature  gradients,  results  in  the  plasma  being 
dragged  up  and  down  the  field  lines.  Hox^’ever,  the  neutral  wind 
is  most  pronounced  at  ionospheric  heights  and  is  not  capable  of 
transferring  ionisation  all  the  way  to  the  protonospherc . It  can 
remove  plasma  to  heights  X\d)erc  tlie  recombination  loss  coefficient 
is  different,  thus  cnhajicing  or  reducing  the  plasma  pressure 
necessary  for  the  diffusion  mechanism. 

3.  Diffusion  of  ))lasma  along  the  i'coiiiagncti  c field  linos  is  the 
most  efficient  way  in  vdiich  plasma  can  travel  up  and  down  in  the 


plasmaspliere  in  rc.sponse  to  pressure  changes.  Ditfnsivo.  motion 


is  not  particularly  important  at  F-layer  heights  because  of 
reactions  caused  by  the  random  movement  of  the  dense  plasma 
but  at  higher  altitudes  the  diffusive  mechanism  is  of  great 
importance  in  plasma  interchange  between  ionosphere  and  protono- 
sphere because  the  ion  gyrof rcquency  exceeds  tlie  collision 
frequency  with  neutral  particles. 

Eetvi7een  600  km  and  1000  km  there  is  a transition  region  at 
which  charge  exchange  occurs  between  oxygen  and  hydrogen.  That 
is,  0 + H <-  0 + H , so  that  the  ion  flux  is  predominantly  0 
below  and  above  this  transition  height.  BANKS  and  DOUPNIK  (1974) 
calculated  these  fluxes  and  found  that,  with  a plasma  density 

— 3 

exceeding  9 x 10^  cm  in  a flux  tube  above  3000  km,  plasma 
could  be  readily  transferred  to  F-region  heights.  Thus  the 


protonosphere  acts  as  a reservoir  of  ions  which  drains  and  refills 
through  interaction  viith  the  ionosphere  at  both  ends  of  the  flux 
tube. 

The  plasmaspheric  flux  has  been  measured  on  a limited  basis 
using  whistlers  (PARK,  1970)  and  incoherent  scatter  (EVANS,  1975a) 
but  continuous  observations  allowing  systematic  study  of  the 
diurnal,  and  seasonal  variations  are  lacKing. 

The  ATS-6  experiment  provides  a method  of  estimating 
integrated  net  protonospheric  fluxes  on  a continuous  basis  as  the 
temporal  change  of  protonospheric  content  is  indicative  of  the 
filling  and  draining  of  this  region  where  production  and  loss 
processes  are  unimportant.  However,  it  should  be  noted  that  the 
protonospheric  content  measured  in  the  ATS-6  experiment  refers  to 
the  slant  path  to  the  satellite  which  intersects  a wide  range  of 
L-value  flux  tubes,  as  shown  in  Fig.  28  for  Aberystwyth  geometry 
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I 

up  to  L -•  4.  The  protonoapher ic  flux  ostiniated  from  AlS-6 
mensurenu-nts  is  thus  an  integrated  parameter  over  tlicse  tubes 
in  contrast  to  the  more  direct  measurements  of  f]ux  essentially 
along  a particular  field  tube  given  by  the  wliistler  and  incoherent 
scatter  techniques. 

b . Results  and  Discussion 

The  monthly  median  diurnal  values  of  have  been  used  to 
estimate  the  temporal  rate  of  change  of  this  parameter  (dNp/dt) , 
calculated  for  each  hour,  which  give's  an  indication  of  the  net 
integrated  ionisation  flux  in  and  out  of  tlie  protonosphere.  To 
reduce  the  scatter  in  the  data  points  the  data  were  first 

filtered  using  a 3-hour  running  mean  before  calculating  the 
gradients.  The  results  for  the  nine  montlis  of  observation  are 
presented  in  Fig.  29  where  positive  values,  indicative  of  an 
upwards  flux  filling  the  tubes,  are  found  for  the  most  part  by 
day  with  negative  values,  corresponding  to  dov^nwards  dei)letion, 
generally  more  prevalent  at  night.  Mention  raust  be  made  of  the 
negative  values  found  during  daytime  in  March.  As  has  already 
been  noted  (Fig.  14),  tlie  median  Np  values  for  these  hours  are 
based  on  only  a small  sample  of  daily  values  at  a time  of  high 
geomagnetic  activity  and  so  must  be  treated  with  caution. 

Similarly  the  impression  of  a strong  dowmwards  flux  during  mid- 
morning in  April  is  in  reality  based  on  only  one  data  point. 

In  an  attempt  to  explain  some  of  the  observed  features  of 
the  diurnal  variations  of  the  integrated  flux,  the  sunrise  and 
sunset  times  at  ionospheric  height  (300  km)  along  the  ray  patli 
estimated  from  the  curves  of  COLIN  et  al.  (1966)  arc  ii\dicated 


h 
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sunset  times  at  300  km  in  the  roiijuj;ali:  ionosplicrc  (36°S,  20*^)^) 
arc  a] so  shown.  Tlie  non-alignment  of  the  geographic  and  geomagnetic 
axes  results  in  sunrise  and  sunset  times  at  the  two  points  conjugate 
to  each  other  which  differ  progressively  for  different  months 
from  a sunrise  coincidence  in  April  and  a sunset  coincidence  in 
February.  It  must  be  remembered  in  interpreting  this  presentation 
that  the  ray  path  crosses  many  field  lines  at  different  heights, 
the  ionospheric  feet  of  which  become  sunlit  at  different  times. 

Fig.  28  gives  an  indication  of  the  range  of  latitude  of  the 
ionospheric  terminations  of  the  field  lines  up  to  L = 4 , but 
taking  account  of  the  longitude  range  as  well,  the  spread  of 
sunrise  and  sunset  times  for  the  feet  of  the  relevant  field  lines 
in  both  local  and  conjugate  hemispheres  can  exceed  ± 1 hour  about 
the  times  shown  in  Fig.  29. 

A second  factor  to  bear  in  mind  is  the  time  constant  between 
an  increase  in  ionospheric  ionisation  and  the  appearance  of 
enhanced  protonospheric  plasma.  A time  constant  'v  1 hour  has  been 
suggested  for  the  diffusive  barrier  at  the  0^/H  transition  height 
(KOHL,  1966),  while  the  travel  time  of  the  ions  to  the  high 

nrotonnsnherp  even  at  thermal  «;r>epd  tnav  rnn<;i  derahl  v exceed  1 hour 

depending  on  the  field  line  path.  Thus  the  overall  time  constant 
may  be  of  the  order  of  several  hours. 

It  might  be  expected  that  changes  in  plasmapausc  position 
would  be  reflected  in  the  observed  data.  For  quiet  magnetic 
conditions  the  plasmapause  has  a minimum  radius  (L  "v  3)  around 
0600  LT  with  a bulge  to  L 'v  5 in  the  evening  sector  (1800  - 2400  LT) 
(CARI’ENTLR,  1966).  The  consequences  to  the  protonospheric  ctuitcnt 
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shape  of  Che  flux  variation,  quickly  inaxiDiising  downwards  v;it)i 


subsequent  i educed  values,  appears  to  indicate  that  the  short 


inner  tubes  are 


Unfortunately  the  curves  for  other  months  do  not  show  such 


clear-cut  trends  and  their  interpretation  is  more  difficult.  The 


downwards  flux  onset  near  either  local  or  conjugate  sunset  is 


apparent  for  all  months,  while  earlier  times  for  November  and 


April  may  be  a reflection  of  the  scatter  of  the  raw  data  points 


It  should  be  noted  that  March,  November  and  April  in  tliat  order 


are  months  of  highest  mean  geomagnetic  activity  in  the  period 


An  upwards  flux  associated  with  sunrise  in  the  local  hemisphere 


can  be  seen  at  all  months  but  it  is 


short-lived  secondary  peak  is  found  around  conjugate  sunrise  in 


December  and  January  wlien  the  time  separation  between  sunrise  in 


the  two  hemispheres  is  a maximum 


It  can  be  concluded  from  these  studies  that  both  local  and 


conjugate  ionospheres  play  a role  in  the  filling  and  depletion 


of  the  protonospheric  flux  tubes  and  that  in  the  interpretation 


of  data  referring  to  observations  along  the  line  of  sight  to  a 


Kr>f  h hr>Tm*  cnhr>Y*r>c 


In  an  attempt  to  gain  further  understand! n; 


interpretation  of  the  protonospheric  observations  the  ionospheric 


f F2  data  for  Slough,  representing  local  conditions,  and  Cai)c. 


Town  (for  a previous  solar  epoch),  appropriate  to  tlie  conjugate 


ionosphere,  arc  plotted  in  Fig.  30.  For  comparison  the  Cape  To\vn 


data  have  been  time  shifted  to  refer  to  the  local  mean  time  at 


Slough.  In  northern  hemisphere  winter  (November  to  February) 


the  f F2  diurnal  variations  at  Slough  show  a single-peaked 


structure  with  maximum  values  of  “v  6 MHz  around  midday.  A similar 


u 


siligl e-pcnlu'd  Viuidtion  is  found  in  the  Cope  Town  data  for  these 
liionths  vnth  waxiiiium  values  n,  7 MMx.  Precise  magnitudes  should  be 
treated  with  caution  as  this  data  refers  to  an  earlier  solar  cycle. 


I 
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However,  the  duration  for  which  f K2  is  in  excess  of  say  5 Ml;z  at 
Cape  Town  is  almost  twice  tliat  at  Slough.  Thus,  from  shortly 
after  local  sunrise  the  protonospheric  content  is  maintained  with 
upwards  fluxes  from  both  hemispheres.  That  t]ie  maximum  fluxes  are 
found  in  general  just  before  noon  is  possibly  a reflection  of  the 
saturation  effect  of  the  lower  L-  shell  tubes  which  have  been 
partially  filled  from  the  conjugate  hemisphere  earlier.  Another 
factor  contributing  to  the  relatively  small  values  of  the  daytime 
fluxes  could  be  the  low  plasma  densities  found  in  the  local 
ionosphere  just  before  dawn  so  that  considerable  time  is  required 
for  sufficient  build-up  of  plasma  pressure  to  produce  significant 
upwards  fluxes.  The  onset  of  dowm;ards  fluxes  in  late  afternoon 
follows  the  droii  in  local  and  to  a lesser  extent  conjugate  f F2 
values  so  that  for  winter  daytime  fluxes  at  least, it  would  appear 
that  the  dominant  interaction  is  that  between  ionosphere  and 
protonosphere.  in  the  local  hemisphere,  llov/ever,  at  night  the 
conjugate  electron  densities  are  somewhat  higlier  than  in  the  local 
ionosphere,  in  December  and  January  at  least,  so  that  a relatively 
strong  upwards  flux  quickly  sets  in  around  conjugate  sunrise.  It 
should  be  noted  that  conjugate  sunrise  will  be  even  earlier  than 
the  time  shciwn  in  Fig.  29  on  other  field  lines  intersecting  the 
path  from  ATS-6  stationed  east  of  Aberystv?yth.  However,  this  net 
upwards  flux  is  short-lived,  the  effects  of  the  predawn  bite  out 
in  fjjl'2  in  the  local  ionosphere  more  t'nan  compensating  so  that 
a net  decrease  in  protonospheric  content  is  found  until  after 
local  sunrise. 


77. 


The  uncorta i I'.Ly  of  some  of  the  March  data  and  the  scatter 
of  the  observed  fluxes  in  April,  a month  of  above?  averaj^e  geo- 
magnetic activity,  mean  Chat  the  fluxes  for  these  months  cannot 
be  interpreted  in  detail.  Consistent  integrated  flux  patterns 
are  found  for  the  northern  hemisphere  summer  months  of  May, 

June  .and  July.  The,  strong  upwards  flux  following  the  earlier 
local  sunrise  quickly  saturates  and  is  follovjed  by  depressed  or 
even  negative  values  around  midday  when  the  local  f^FZ  values 
are  kept  almost  constant  under  the  dynamic  influence  of  the 
neutral  wind.  The  evening  peak  in  the  local  ionosphere  is  more 
than  compensated  by  the  rapid  drop  in  electron  density  in  the 
short  winter  day  at  the  conjugate  point  so  that  strong  downwards 
fluxes,  it  would  seem  to  the  conjugate  ionosphere , arc  observed. 

These  quickly  maximise  then  drop  indicating  an  emptying  of  the 
inner  L-  shell  tubes. 

It  is  thus  apparent  that  an  interpretation  of  the  observed 
protonospher ic  variations  requires  consideration  of  interactions 
with  the  ionosphere  in  both  hemispheres  and  while  in  the  main  the 
integrated  protonosphcric  content  is  determined  by  local  conditions. 


there  arc  times  when  the  conjugate  ionosphere  is  the  dominant 
factor . 

1 

1 Average  integrated  fluxes,  presented  in  Table.  7 for  each 

month,  show  that  to  within  the  accuracy  of  the  estimates  the 
magnitude  and  duration  of  the  upwards  flux  by  day  is  broadly 
compensated  by  the  downwards  flux  at  night.  PARJC  (1970)  has 

I calculated  average  daytime  and  nighttime  fluxes  from  a ]imit;ed 


data  base  of  whistler  observations  of  tube  content.  He  estimates 
a flux  through  the  1000  km  level  with  a daytime  filling  magnitude 


Table  7.  AVEKAGK  i:Tl'ECrb\Tl,l)  TROTOEOSTlIi.iU  C I'LDXEb 


Upwards 

Downwards  j 

7 _?  -1 

X 10'  era  r> 

Duration  hr 

X 10^  cm  s 

Duration  hr 

N 

Nov 

2.2 

12 

-2.2 

12 
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2.3 

12 

-3.0 

11  f 

I Jan 

i.9 

lA 

-2.8 

10 

Feb 

2.9 

11 

-2.5 

13 

Mar 

4.0 

7 

-4.3 

11 

Apr 

5.4 

10 

-4.2 

13 

May 

4.2 

15 

-6.1 

9 

1 

1 Jun 

3.8 

9 

-4.2 

10 

Jul 

3.9 

8 

-3.4 

13 
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factor  of  two.  These  values  are  si.giuf  i cantly  higlicr  th.m  llie 
results  of  the  present  study.  However,  I’ARK's  data  refer  to  a 
quiet  period  du)ing  which  the  protonosphore  drains  and  refills  in 


a regular  diurnal  pattern.  Secondly,  his  results  correspond  to 
1000  km  level  while  the  ATS-6  observations  are  for  a flow  at 
about  2000  to  lOOO  km  where  the  fluxes  are  probably  smaller  and 
the  cube  area  is  increased  by  a,  5X.  Another  important  factor 
is  that  the  observations  of  PARK  refer  to  tubes  for  L > 4 while 
the  Abery-stv.^th  ray  path  crosses  tubes  for  L ^ 1.7.  CARPENTER  and 
PARK  (1973)  hi’ve  argued  that  the  flux  is  enhanced  for  higlier  L- 
sliolls  because  of  Che  relatively  less  dense  plasma  in  these 
larger  volume  tul>es.  Finally,  as  will  be  seen  later,  the  protono- 
spheric  content  builds  up  consistently  during  a long  quiet  period 
with  the  daytime  flux  exceeding  that  at  night.  This  trend  is 
interrupted  sharply  during  moderate  to  severe  geomagnetic  activity 
and  many  days  or  even  week.s  are  required  for  recovery.  Thus  on 
average  the  daytime  and  nighttime  fluxes  do  not  differ  significantly 
in  magnitude. 

8.  SL,\NT  Sl,AF>  THICKNESS 


Slab  thickness  is  a parameter  which  give.s  an  indication  of 
the  shape  of  the  ionosplieric  profile  by  normalising  the  electron 
content  with  respect  to  the  peak  electron  density  • Following 

the  notation  of  DAVIES  et  al.  (1976)  the  slant  slab  thickne.ss 
(t^)  is  given  by 


t 

s 


1 m 


16. 


Values  of  can  be  obtained  from  published  r2-  layer  critical 
frequency  data.  For  comparison  of  .slant  slal)  thicknc.s.s  with  the 


otlier  measuri  d fdiapc  parameter,  the  F-  f. actor,  it  is  necessary  to 
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evaJuaLt'  at  the  sub-ionospheric  pf)iiiL  of  the  ray  path  ( 47.'}" 

3.3"h  for  35->  hm  altitude).  However,  a.s  there  i .s  no  convtuii  ent 
ioiio.sonde  at  this  location  an  extrapolation  procedure  has  beer, 
used  to  generate  the  required  data  from  values  obtained  from  a 
number  of  ionosondes  in  the.  general  locality. 

A prior  tost  of  fyF2  data  from  several  European  stations 
showed  a substantial  longitudinal  gradient  in  the  values  for 
the  same,  local  mean  time.  Therefore,  only  four  stations,  with 
the  smallest  range  in  longitudes  about  the  sub-ionospheric  point, 
have  been  used:  De  Bilt  (.52.i"n,  5.2"e),  Slough  (SI.S'^'n,  0.6"w)» 
Dourbes  (50.i"n,  4.6”e)  and  Lannion  (48.7"n,  3.5°W).  Fig.  31 
show,  the  geographic  locations  of  these  stations  with  respect  to 
the  ground  proicction  of  the  ray  path.  Data  w'as  unavailable 
for  Lannion  for  March  1976  and  for  De.  Bilt  from  March  1976.  The 
missing  values  luive  been  interpolated  from  the  published  v.vrldv;ide 
contours  of  f F2  (CCIR,  1967)  bv  usinc  the  existino  data  for  the 
other  stations. 

The  monthly  median  fQF2  values  for  each  hour  from  the  four 

stations  have  been  considered  as  a function  of  latitude  and  f F2 

o 

at  the  sub-ionosphcric  latitude  obtained  from  a legression  line 
analysis.  When  tlie  correlation  coefficient  was  less  than  0.7 
the  fQF2  value  at  the  sub-ionospheric  point  was  estimated  from 
the  CCIR  plots.  A comparison  of  the  tw’o  methods  shov;ed  that 
differences  were  in  general  within  the  errors  in  scaling  the 
contour  plots. 

It  should  be  noted  that  to  keep  the  workload  reasonable  in 

this  preliminary  study  attention  was  concc-ntratcd  on  monthly 

median  f F2  and  N„,  values  so  that  the  'median'  t:)ah  thickness 
o 1 


calculated  may  differ  slightly  from  the  true  median  of  tlie  individual 
daily  values,  however,  it  is  unlikely  t!iaV  any  conclusions  are. 


affected  by  this  approximation.  In  fact,  following  DAVTliy  et  al. 
(1976)  a further  smoothing  has  been  carried  out  on  the  data  by 
calculating  a 3-  hour  running  mean  of  the  hourly  values  for  each 
month . 

The  smoothed  median  values  of  slant  slab  thickness  are 
presented  in  Fig.  32.  Daytime  values  < ^00  km  in  winter  give  way 
to  thicknesses  'v  800  km  in  summer  while  by  contrast  magnitudes 
well  in  excess  of  700  km  on  winter  nights  contrast  to  thicknesses 
< 600  km  during  darkness  in  summer.  Also  of  note  ore  the  predawn 
and  post  sunset  peaks  in  slab  thickness  which  are  most  prominant 
in  the  winter  months. 

Following  DAVIES  ct  al.  plots  have  been  made  of  F against 
for  each  month  (Fig.  33).  kTiilc  the  daytime  values  of  F are 


the  predawnn  sector  increases  in  F are  matched  by  l.orge  slab 
thickness  in  v;inter  when  tlie  protonospheric  content  assumes 
greatest  proportional  importance.  The  curves  show  essential 
similarities  to  those  of  DAVIES  et  al.  being  'open  loops'  in 
sunmier  and  more  nearly  linear  in  winter,  although  the  secondary 
loop  found  near  midnight  in  winter  reflects  the  decrease  in 
slab  thickness  observed  in  the  European  data  and  not  in  the 
American  values. 


9 . COMPARISON'  ITil  AlCiRlCAN  .SECTOR  DATA 
a.  Inti'oduction 

Although  the  prime  ol-.jcctivc  of  the  present  v:ork  has  been 
the  study  of  i onospheri c/protonospher ic  coujil  ing  along  the  .ATS-6 
ray  path  to  Aberystwyth  it  was  considered  important  that  the.  data 


F'i. 


BDSSSCSSESSBS?^ 


Slant  Slab  Thickness 


obtaiiicul  sliuuld  be  tonipared  to  other  available  ATS-0  results. 


FRI'I’/’,  (1976)  lias  [uibl  Lshed  a complete  set  of  data  for  the  Boulder, 
USA  station  (4o‘^’n,  lOS^^W)  obtained  from  duly  1974  to  May  1975 
during  the  first  phase  of  the  ATS~6  project  when  the  satellite 
was  stati  ned  at  94°W  longitude.  Apart  from  individual  case 
studies  of  a few  days  duration  and  theoretical  model  work  this 
experimental  data  forms  the  only  base,  at  the  time  of  writing, 
with  Vv’hich  to  compare  the  Aberystviyth  results.  Surprisingly, 
at  first  sight  at  least,  the  monthly  median  Boulder  protonospheri c 
content  shows  a diurnal  variation  in  which  minimum  values  are 
found  by  day  with  a maximum  at  night,  especially  for  the  winter 
months.  This  pattern  in  American  sector  data  has  been  repeated 
in  the  case  studies  of  SOICHER  (1976b)  and  P0LETTI-LIL71ZI  (1976), 
but  has  been  passed  over  largely  without  comment  by  all  three 
groups  of  Workers.  The  very  limited  data  to  come  as  yet  fi'om 
the  European  sector  vjovld  appear  to  be  in  essential  agreement 
with  the  Aberystwyth  results  in  terms  of  I he  form  of  the  diurnal 
variation  (SOICHliR,  1976b  and  HARTMANN’  et  al.  1976). 

It  thus  appears  that  there  is  a discrepancy  betvjecn  American  and 
European  observations  of  protonospheric  content  which  requires 
expl anat ion . 

b.  Comparison  of  Results 

Tlic  publislicd  Boulder  monthly  median  data  for  the  ATS-6 
Phase  I period,  converted  to  local  mean  time,  are  compared  to 
correspond iitg  monthly  values  for  Aberystv.ylh  from  the  later 
Phase  ir  ope.ratii'us  in  Fig.  34.  The  differences  in  the  form  of 
the  diurn.'’]  variations  and  in  the  N.^,  magnitudes  for  the  two 


stations  aro  immediately  apparent.  An  examination  has  shown  tliat 

those  differences  cannot  be  attributed  to  differences  in  the  level 

of  solar  and  geomagnetic  activity  between  the  data  sets.  However, 

a table  corresponding  to  Table  1 for  Aberystvjy th,  giving  the 

geophysical  parameters  for  the  ATS-6  to  Boulder  geometry  allows 

certain  conclusions  to  be  drawn.  Firstly,  the  discrepancies  in 

Np  cannot  be  attributed  to  range  differences  along  the  slant  path. 

Wliile  the  range  through  the  ionospheric  part  of  the  path  up  to  a 

height  of  2000  km  is  m 2600  km  for  Boulder  compared  to  m 3600  km 

for  Aberystvjyth,  the  range  difference  of  the  protonospheric  paths 

is  only  m 3%.  Secondly,  for  the  Boulder  geometry  fj  is  a single 

Fig. 35:  . 

valued  function  of  height/  the  height  increase  in  cos  0 is  too 

small  to  offset  the  decrease  in  gyrof requc.ncy  as  is  the  case  for 

the  Aberystu’yth  geometry.  As  a result,  f^  for  Boulder  is  a much 

more  sensitive  function  of  height  around  400  km  than  that  for  . 

Aberystwyth,  so  that  a change  of  + 100  km  in  the  height  at  which 

f.  is  evaluated  can  result  in  a change  of  + 50%  in  K., . Hence 

in  the  interpretation  of  the  Boulder  diurnal  variation  of  Np  it 

is  of  great  importance  to  assess  the  effect  of  di.urnal  changes 

in  the  height  of  the  F2-  layer.  A procedure  similar  to  that 

already  described  for  the  Aberystwyth  data  was  carried  out  to 

estimate  the  true  height  of  the  ionospheric  peak  (h^)  from 

ionosonde  data.  Unfortunately,  the  World  Data  Centre  were  unable 

to  provide  f F2  and  f E values  for  either  Boulder  or  WTii te  Sands 

for  the  period  of  Phase  I ATS-6  observations.  However,  monthly 

median  ionosonde  data  from  Boulder  for  corresponding  months  of 

the  jirevious  solar  cycle,  eleven  years  earlier,  have  been  used  to 

estimate  h and  its  diurnal  variation, 
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The  substanLinl  effect  of  this  varintion  of  h on  tlic  values 

in 

of  particularly  for  daytime,  is  il  1 ustr/iteci  in  Fij;;;.  36  - 38. 

Clioosine  h + 100km  as  the  heieh-t  for  the  estimation  of  f,  it  can 
m L 

be  seen  from  Fig.  39  that  the  afternoon  values  of  are 
significantly  enhanced.  For  the  equinoxial  month  of  March  a niinimiim 
value  of  Np  is  now  found  around  dawn  rather  than  the  broad  morning 
minimum  of  the  raw  data  for  which  a constant  f evaluated  at  about 

[j 

380  km  had  been  used.  In  winter  a raising  of  the  daytime  values 
is  detected,  but  this  still  leaves  them  lower  than  those  found 
at  night  and  in  fact  for  no  reasonable  choice  of  increment  to  li 
can  the  diurnal  variation  be  reversed.  Thus,  while  the  remarks 
of  KLOBUCHAR  and  KIDD  (1972)  and  ALMDIDA  (1974)  are  partially 
confirmed  and  the  diurnal  variation  of  the  layer  height  plays  a 
m.ore  significant  role  for  stations  in  the  American  sector  tlian 
Europe,  nevertheless,  an  alternative  ionospheric  as  opposed  to 
observational  explanation  must  be  sought  to  account  for  the  witiLei 
daytime  minimum  in  Np. 

Differences  in  the  magnitude  of  Np  for  Abcryst\gyth  and 
Boulder  also  require  interpretation.  It  can  be  seen  from  Fig.  34 
that  while  nighttime  values  in  vjintcr  are.  comparable,  at  other 
hours  and  seasons,  in  particular  sumaser,  the  Boulder  content 
may  be  less  than  half  that  found  at  Aberystwyth. 

Consideration  of  the  ray  path  geometry  with  respect  to  the 
field  lines  for  Boulder  (Fig.  40)  shows  that  at  comparable  points 
on  the  path  the  intersection  is  with  higher  L-sliell  tubes,  with 
a minimum  I.  'v  2,  than  for  Aberystwyth  (Fig.  28)  where  a 
considerable  [Portion  of  the  path  interseets  field  lines  for  which 
L < 2.  In  addition  the  field  line  geometry  result;;  in  a con,i<igaLe 
ionospheric  point  at  much  higher  latitude  than  Boulder  itself  in 
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contrast,  to  the  lower  conjuj’.ate  latitude  for  Aberystv^^th.  Thus, 

(UKLLCHAM'IKRS  and  PIGGOTT,  1960), 
remembering  the  anomalous  behaviour  of  the  Antarctic  ionosphere/ 

it  seems  possible  that  the  explanation  invoked  to  account  for  the 

Aberystvjytli  results,  involving  i onospheric/protonospheric  interaction 

at  both  ends  of  a field  line,  may  provide  some  qualitative 

interpretation  of,  not  only  the  very  low  values  in  summer  but 

also,  the  inverted  diurnal  variation  in  the  Boulder  data  for 

wi  liter. 

To  examine  this  further,  the  monthly  median  fjjF2  data  for 
Boulder  mentioned  earlier  have  been  plotted  in  Fig.Al  together 
with  data  from  a similar  period  from  Argentine  Island  (65*^S,  64*^W) , 
this  station  being  the  ionosonde  closest  to  the  location  (57°S, 

127°W)  conjugate  at  300  km  to  tlie  Boulder  sub-ionospheric  point. 

The  latter  data  have  been  time  sliifted  to  correspond  to  tlie 
local  mean  time  at  Boulder. 

As  an  additional  aid  to  interpretation,  the  integrated  fluxes 
(dXp/dt)  have  been  computed  for  Boulder  using  3-hour  running  mean 
smoothed  Nj,  data  corrected  by  use  of  a variable  f^  evaluated  for 
h^  + 100  km.  These  are  plotted  in  Fig.  42. 

Considering  first  July,  upwards  fluxes  are  observed  to 

commence  soma  2 to  3 hours  after  local  sunrise  when  f F2  for 

o 

Boulder  rises  from  a minimum  value  about  3Mllz  to  'v  SMHs  before 
noon.  The  integrated  fluxes  are  smaller  than  for  Aberystwyth  at 
this  time  and  do  not  shov;  quite  the  same  rapid  saturation  because 
the  average  L-  value  of  the  intersecting  tubes  is  greater.  The 
small  magnitude  of  the  upwards  fluxes  and  the  low  Np  values  for 
July  can  be  explained  in  terms  of  the  long  winter  night  of  the 
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relatively  liigh  latitude  conjugate  ionosphere,  with  for 

Argentine  Island  < 2.51111?.,  for  some  20  hours.  Thus  in  northern 

hemisphere  summer  there  is  a steady  drain  of  plasma  to  the 

conjugate  ionosphere  for  most  of  the  day  with  little  upwards 

compensation  even  during  the  short  daylight  hours.  It  is 

interesting  to  note  that  in  July  the  £^F2  diurnal  variations 

for  Slough  and  Boulder  are  very  similar  in  form  and  magnitude 

but  the  markedly  different  Np  data  for  the  tv70  stations  can  be 

accounted  for,  at  least  qualitatively,  by  the  differences  in 

conjugate  ionospheric  behaviour  and  path/flux  tube  geometry. 

For  December,  the  noteworthy  feature  is  the  anomalous 

behaviour  of  fQF2  for  Argentine  Island  in  response  to  the  300  km 

ionosphere  being  sunlit  24  hours  a day.  In  fact,  it  should  be 

noted  that  even  at  the  Boulder  conjugate  latitude  (57°S)  there 

is  no  ionospheric  sunset  from  November  to  January.  The  resulting 

f F2  values  minimise  at  'v  5 MHz  in  early  afternoon,  under  the 
o 

influence  of  the  polewards  neutral  v.’ind,  (KOHL  and  KING,  1967), 
but  by  night  the  F2  critical  frequencies  become  "v  8MHz. 

The  Boulder  f^F2  values  range  from  a,  3 MHz  at  night  to  7 Ifllz 
at  the  single  peak  maximum  in  the  middle  of  the  winter  day.  The 
much  higher  Np  values  for  December  compared  to  July  can  be 
attributed  to  the  generally  higher  electron  densities,  particularly 
in  the  conjugate  ionosphere.  In  addition,  the  daytime  miniwum  in 

Np  reflects  the  anomalous  behaviour  of  tlie  conjugate  iono.sphere 

whose  influence  is  dominant  in  northern  hemisphere  winter  daytime. 

It  can  be  noted  witli  caution  that  use  of  h + 150  km  in  the 

m 

correction  of  the  Np  data  raises  tlie  winter  midday  values  further 
removing  the  minimum,  a situation  wiiich  may  take  account  more 
closely  of  the  behaviour  of  the  local  ionosphere. 
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IL  is  possible  to  consider  e.ich  iiiontli  in  detail  and  interpret 


i 


features  of  both  the  ilata  and  the.  intcj^ratcd  pet  fluxes  in 
tcrnis  of  ionospheric,  behaviour  in  the.  two  hemispheres.  However, 
wc  shall  discuss  only  March  since  for  this  month  the.  conjugate 
ionosonde  data  sliov/s  a dramatic  change  to  a daytime  peak  v.dth  the 
arrival  of  ionospheric  darkness  in  March,  but  the  data  after 
correction  h.tve  a dawn  minimum  and  a nighttime  peak  which  at  first 
sight  at  least  appears  spuriou.s.  Some  2 to  3 hours  after  sunrise 
small  positive  fluxes  are  found  which  persist  (apart  from  one 
possibly  anomalous  point  although  even  this  may  be  interpreted 
by  a more  complete  model  of  the  actual  physical  situation)  to 
well  after  local  and  even  conjugate  sunrise.  Tlie  post  midnight 
strong  upwards  fluxes  and  the  early  morning  peak  in  Np  may 
represent  a protonospheri c filling  from  the  late  evening  peak 
in  the  conjugate  hemisphere,  delayed  by  several  hours 

because  of  the  relatively  long  flux  tubes. 

While  it  is  recognised  that  there  are  limitations  to  this 
fairly  simple  qualitative  approach  to  the  interpretation  of  the 
protonospheri c data,  nevertheless  the  results  obtained  for 
Aberystw^'th,  together  with  some  confirmation  from  the  contrasting 
Boulder  situation,  would  appear  to  justify  further  model  studies 
of  a quantitative  nature.  These  should  take  account  of  both 
local  and  conjugate  ionospheres  and  use  multiple  flux  tube 
geometry  with  appropriate  time  constants  v.’liich  would  more  closely 
approximate  to  the  actual  experimental  conditions. 


( 
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10.  CASl']  STDDIKS 


a.  Introduction 


This  report  has  been  concerned  so  far  with  monthly  median 
values  of  various  parameters  determined  by  the  ATS-6  radio  beacon 
experiment.  While  these  have  provided  us  with  information  on  tlie 
average  beliaviour  of  ionosphere  and  protonosphere  they  give  little 
indication  of  the  response  to  individual  events,  in  particular  to 
changing  levels  of  geomagnetic  activity.  The  response  of  ionospheric 
electron  content  to  geomagnetic  storms  is  reasonably  well  documented 
(e.g.  MENDILLO  and  KLOKUCHAR,  1974)  so  atte.ntion  has  been 
concentrate.d  here  on  the  day  to  day  behaviour  of  protonospheric 
content.  As  a short  conclusion  to  this  preliminary  study  of  ATS-6 
data  two  periods  have  been  selected  for  discussion,  both  cover 
t?  f major  storm  activity  but  one  also  takes  in  a two  week 


liet  days  during  which  rarely  exceeds  3. 


b.  13  November  - 4 December,  1975 


Hourly  values  of  Np  available  throughout  this  period  are 

plotted  in  Fig.  43,  Also  showna  for  comparison,  recurring  for 

each  day,  is  a plot  of  the  average  Np  for  a 30  day  period  centred 

on  22  November  1976.  Values  of  the  3-hour  K index  are  presented 

P 

above  the  corresponding  Np  data. 

The  period  starts  with  four  days  vjhen  Kp  never  exceeds  3 and 
the  Np  values  are  seen  to  rise  gradually.  Enhanced  geomagnetic 
activity,  maximising  on  17  November  causes  a drop  in  Np  around 
1800  LMT  on  tliat  day  and  values  remain  low  during  tlie.  moderate 
activity  period  to  20  November.  Following  = 5_  early  on 
21  November,  Np  maximises  around  midday  and  relatively  high  values 
are  maintained  on  22  November.  A sudden  commencement  at  2306  UT 
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on  21  November,  witli  K = (> . on  the  uf  teriiocni  of  the  next'  dav,  is 

p ■) 

folJov.'ed  by  a sharp  drop  in  Kj,  around  1800  UT  on  22  Novinnber.  A 
steady  refilling,  of  Llie  protonospiiero  is  found  till  the  29  November 
with  a clear  diurnal  pattern  to  be  seen  in  N^,.  Kidianced  values  of 

Kp , reaching  5 on  29  November  are  accompanied  by  high  N^  on  that  day 
followed  by  a gradual  dcjjlct ion  superposed  on  the  regular  diurnal 
variation  till  the  morning  of  2 IJecember  while  remains  in  the 
main  > 4.  The  jirotonospher ic  content  shows  signs  of  subsequent 
recovery  on  3 December  as  geomagnetic  activity  lessens. 

c.  28  April  - 21  May,  ]_97f) 

The  period  (I’ig.  44)  starts  with  four  days  in  which  the  geomagnetic 
activity  reaches  Kp  = 4^  late  on  29  April  and  Np  shows  a fluctuating 
pattern,  particularly  on  30  April  and  1 May  when  the  diurnal  variation 
is  almost  masked  by  the.  apparently  random  variations.  A maior  storm 
with  sudden  commencement  at  1828  UT  on  2 May  and  Kp  reaching  8^  is 
preccided  by  a sharp  increase  in  Np  around  midday  to  be  followed  by 
a depletion  with  Np  reaching  very  low  values  on  the  morning  of  3 May 
and  even  lower  magnitudes  around  midnight  of  that  day.  In  fact,  this 
value  1.3  x 10’®  m war.  tlie  lowest  Np  found  during  the  entire 
observing  period.  Three  days  follow  \cith  substantial  breaks  in  the 
Np  data,  but  the.  values  plotted  shov;  the  trend  of  a replcnisliing 
protonospher ic  content  which  is  continued  until  about  16  May,  during 
a time  of  quiet  geomagnetic  activity.  Superposed  on  this  trend  is 
a regular  diurnal  variation  witli  increasing  amplitude  as  the  quiet 
period  continues.  From  16  to  19  May  geomagnetic  ai’.tivity  is  very 
small  and  Np  is  consistently  above  average  with  a continuing  clear 


larjjo  diurnal  variafion.  IL  would  thus  appear  that.  t!io  deiileted 
j)rolonosphere  of  1 May  requires  some  lA  days  of  quiet  eondit ions 
with  steady  diurnal  partial  filling  and  draining  to  reach  a quasi- 
equilibrium situation  about  16  May.  Thus  from  16  to  18  May  the 
protonospheri c content  for  ATS-6  to  Aberystwyth  geometry  presents 

its  equilibrium  pattern  for  quiet  geomagnetic  activity  with  values 

• • 1 6 

maximising  around  7 x 10  m during  daytime  filling  and 

1 6 

depleting  to  about  4 x 10  m at  night. 

This  pattern  is  confirmed  from  the  presentation  of  Fig.  45 
where  the  maximum  and  minimum  values  of  Np  for  each  day  throughout 
this  period  are  plotted  and  show  that  while  both  increase,  steadily 
from  the  storm  depletion  the  daytime  filling  more  than  compensates 
the  nighttime  draining  for  seme  12  to  14  days.  Further  confirmation 
is  found  in  the  increase  in  daytime  Np  maximum /N^X  where  in  the 
equilibrium  situation  the  daytime  protonospheri c content  contributes 
some  20%  to  the  total  as  opposed  to  only  11%  on  the  day  following 
the  sudden  commencement. 

d.  Discussion 

The  first  point  to  note  from  this  study  of  Np  changes  during 
two  extended  periods  is  that  the  regular  diurnal  variation,  with 
early  morning  minimum  and  afternoon  maximum,  is  essentially  a 
quiet  time  phenomenon.  UTien  values  exceeding  3 are  found  over 
a period  of  several  days  the  diurnal  pattern  of  Np  becor.ies  irregular 
with  noiseliUe  quasi-random  fluctuations  which  tend  to  mask  the 
diurnal  component  of  the  variation. 

Enhanced  geomagnetic  activity  accompanying  a sudden 
commencement  is  associated  with  an  increase  in  the  daytime 
protonosphcric  content.  This  re.sult  is  consistent  with  that 
reported  by  I’OhETTl -hlUZZI  et  al.(1976)  and  agrees  with  the 
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nic.isiireiiK'nUs  of  density  above  '4OO  km  iliscussed  by  IlAI’.b’l.ASHI 
ot  al.  (1976)  who  report  an  ujiwards  flow  of  1 ip, lit  ions  at.  hi>;li 
latitudes  associated  with  the  storm  sudden  corriiencenient . It  is 
well  known  that  ionospheric  electron  content  can  show  a signi- 
ficant positive  effect  during  daytime  following  a sudden 
coimnencement  (KKRSl.EY  and  DAS  GUPTA,  1976).  This  daytime  increase 
is  associated  with  a raising  of  the  ionisation  under  the  influence 
of  the  niagnetosphcric  electric  field  penetrating  to  the  mid- 
latitude ionosphere  and  also  neutral  winds,  possibly  reversed 
in  direction  and  blowing  equatorwar-ds  in  response  to  enhanced 
temperatures  in  the  auroral  zone.  It  is  thus  to  be  expected 
that  the  increase  in  ionospheric  plasma  will  he  reflected  in 
enhanced  protonospheric  content.  7\n  additional  factor  which  may 
cause  some  enhancem.ent  in  slant  is  the  expansion  of  the 
plasmaspherc  in  the  afternoon  sector  in  response  to  increased 
magnctospheric  convection. 

The  rapid  depletion  in  Np  as  tlie  line  of  sight  moves  into 
I the  nightside  is  consistent  with  the  corresponding  drop  in 

ionospheric  electron  content  and  peak  density  reported  by  many 
workers.  The  main  cause  of  this  decrease  is  almost  certainly 
linked  to  changes  in  thermospheric  circulation  reflected  in  an 

'1 

increased  N2/O  ratio  (PROLSS  and  VOX  ZAUN,  1974),  but  the 
enhanced  magnctospheric  electric  field  probably  plays  a role  in 
depressing  the  nighttime  ionisation  and  thus  reducing  the. 
ionospheric  electron  content  (Kl^RSLEY  and  DAS  GUPTA,  1976)  ainl 

r 

hence  protonospheric  content.  An  additional  facLir  vdiich  will 

* 

contribute  to  the  depression  of  nighttime  Np  during  tin;  main 


\ ' 


found  that  tiio  i)lasniapau!.t;  moved  in  to  I,  ==  2.4  .several  Iiourn  after 
local  midnight  during  a storm  in  vdiich  Kp  reached  7. 

The  slow  r(?Fillin;;  ol  the  protono.sphorc  alter  the  storm 
depiction  was  found  to  taVie  14  days  to  reach  saturation  during 
the  quiet  conditions  in  May.  The  results  of  the  whistler  studies 
rej)orted  by  TARK  (1974)  indicate  time  constants  to  reach 
equilibrium  ranging  from  1 day  at  L - 2.5  to  'v  8 days  at  L = 4. 

The  storm  of  2 May  1976  v;ith  Kp  maximising  at  8^  was  more  intense 
than  that  discussed  by  PARK  so  that  the  protonospherie  depiction 
may  have  been  more  complete  in  our  observations  where  the  lowest 
Np  value  ever  found  occurred,  however  in  terms  of  season  and  solar 
minimum  epoch  the  data  arc  comparable.  It  is  interesting  to  note 
that  park's  observations  were  interrupted  by  subsLorm  activity 
after  8 days  so  that  he  was  unable  to  study  possible  saturation 
beyond  1,  “ 4.  However,  in  the  present  observations  the  period 
of  quiet  extends  to  'v  17  days  allowing  complete  equilibrium  along 
the  slant  path  to  be  reached  after  s.  14  days  when  the  plasinapause 
radius  had  extended  to  its  furthest  possible  quiet  time  position. 

An  adtlitional  point  of  interest  is  that  after  v-  7 - 8 days  the 
Np  magnitudes  and  their  diurnal  variation  become  comparable  to 
the  30-  day  average  value.s  suggesting  that  the  median  protonospherie 
content  and  flux  data  discus,sed  earlier  in  tins  report  represent 
the  integration  over  flux  tubes  up  to  L s,  4.  This  observation 
can  be  linked  to  PARK's  comment  that  the  pla5:ma3phere  usually 
con.sists  of  two  distinct  regions;  an  ' nnor'  pi  asmaspberc , which 
is  in  equilibrium  witli  the  underlying  isphcrc  over  a 24-  hour 

period,  and  an  'outer'  jil  asmaspberc  I.  > 4. .5,  which  is  still 

recovering  from  the  previous  disturb;  e. 


w 




The  average  daytime  uiiwnrds  integrated  prot onosi*her i c flux 

for  16  May,  when  equilibrium  conditions  had  been  reached,  was 
^ -2  _ 1 

8 X 10  cm  s which  can  be  compared  to  a filling  llnx 

_ - 2 _ J 

“v  3 X 10  cm  s at  1000  km  for  L > 4 reported  by  PAIIK  (1970) . 

It  should  be  noted  that  PARK  has  Interpreted  his  data  in  terms 
of  interaction  with  the  ionospliere  at  one  end  of  the  flux  tube 
while  the  present  work  suggests  that  the  conjugate  hemisphere 
must  also  be  taken  into  consideration.  Thus  the  present  integrated 
value,  only  some  25%  of  PARK's  would  appear  to  give  added 
confirmation  to  the  suggestion  made  earlier  in  this  report  that 
saturation  of  the  short  inner  tubes  with  L < 2 occurs  after  only 
a few  hours  filling  almost  every  day. 

11.  CONCLUSIONS 

The  work  reported  here  represents  a first  stage  analysis 
of  the  ATS-6  radio  beacon  data  for  Aberystwyth.  Attention  has 
been  concentrated  in  the  main  on  hourly  values  of  \'arious  parameters 
treated  as  monthly  medians. 

A consideration  of  the  overall  accuracy  of  the  modulation 
phase  and  Faraday  rotation  measurements  has  indicated  that  in 
general  the  observational  errors  result  in  absolute  protonospher ic 
content  data  accurate  to  about  ± 15%.  The  short  term  relative 
errors  are  usually  much  smaller  so  that  a reliable  picture  of 
dynamic  protonospher ic  changes  can  be  obtained  from  ATS-6 
observations.  It  has  been  shown  that  for  geometry  applicable 
to  the  Aberystwyth  measurements  the  longitudinal  gyrof requency 
used  in  Faraday  content  evaluation  is  relatively  insensitive 
to  height,  so  that  the  error  in  resulting  from  neglect  of  tlie 
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IS  in  gene 


wns  found  to  be  i.iucb  more  favourable  than  for  a station  in  the 


American  sector  where  the  corresponding  error  could  be  ns  large 


Tbus  correction  for  layer  height  changes  is  of  critical 


importance  in  the  analysis  and  interpretation  of  N data  for  tlie 


American  stations  observing  during  Phase  I ATS-6  operations  but 


much  less  so  for  inea 


urements  made  in  Europe  during  Phase  II 


The  dynamic  interchange  of  ionisation  be-tween  ionosphere  and 


protouospbore  has  been  studied  by  estimating  hourly  values  of 


average  integrated  fluxes  from  the  temporal  gradient  of  protono 


spheric  content.  Tlie  protonospbcric  flux  shows  a coiiiplex  behaviour 


both  diurnally  and  seasonally,  but  some  considerable  progress  has 


been  made  in  the  interpretation  of  the  results  on  a qualitative 


basis  by  considering  diffusive  interaction  between  ionosphere 


and  protonosphere  in  both  local  and  conjugate  hemispheres.  This 


interpretation  has  been  given  added  confirmation  by  application 


to  totally  contrasting  protonospbcric  data  from  an  American 


station  which  nevertheless  can  he  related  to  ionospheric  behaviour 


in  the  tv:o  hemispheres 


The  Aberystwyth  observations  indicate  that  the  inner  flux 


and  that  the  average  behaviour  of  N„  reflects 


hours  almost  e 


a diurnal  partial  draini.ng  Vjut  complete  refilling  of  tube 


which  L < 4.  Following  a geomagnetic  storm  depletion  some  7 days 


are  required  to  saturate  these  tubes,  but  approximately  14  quiet 


days  are  necessary  to  reach  an  absolute  dynamic  satur.ition  along 
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